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ABSTRACT 

We present forty-one near infrared (NIR, 0.7-2.5 fim) spectra from normal 
Type la supernovae (SNe la) obtained at epochs ranging from fourteen days 
before to seventy-five days with respect to the maximum light date in the V- 
band. All data were obtained at the IRTF using the SpeX instrument. We 
identify many spectral features, measure the Doppler velocities, and discuss the 
chemical distribution of explosion products in SNe la. We describe procedures 
for smoothing data, fitting continua, and measuring absorption features to insure 
consistency for measurement and analysis. 

This sample provides the first opportunity to examine and compare a large 
number of SNe la in this wavelength region. NIR data are a rich source of 
information about explosion products whose signatures are blended or obscured 
in other spectral regions and NIR observations probe a greater radial depth than 
optical wavelengths. We analyze similarities and differences in the spectra and 
we show that the progressive development of spectral features for normal SNe la 
in the NIR is consistent with time. We confirm the presence of O I, Mg II, Ca II, 
Si II, Fe II, and Co II in these SNe. Possible identifications are made for S I, 
Si III, Mn II, and Fe III. There is no evidence in these data for H I, He I, C I, or 
CII. 

As the explosion products expand and cool, progressively deeper layers are 
revealed. Thus a time sequence of spectra examines the chemical structure and 
provides direct evidence of the physical properties of SNe la from the outer 
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layers to deep inside the SN. Measured Doppler velocities indicate that burning 
products in SNe la are distributed in distinct layers with no large scale mixing. 
Carbon is not detected in these data, in agreement with previous results with 
NIR data establishing very low limits on carbon abundance in SNe la. Carbon 
burning products, O and Mg, are plentiful in the outer layers suggesting that 
the entire progenitor is burned in the explosion. The data provide a resource 
for investigations of cross-correlations with other data libraries that may further 
constrain SN la physics and improve the effectiveness of SNe la as cosmological 
distance indicators. 



Introduction 



This sample of forty-one near infrared (NIR, 0.7-2.5 /xm) spectra from Type la super- 
novae (SNe la) provides the first opportunity to study NIR spectral evolution in a large 
number of SNe la. 

SNe la are intrinsically interesting because they are very powerful explosions (~ 10^^ 
ergs) involving extreme and exotic physics. SNe la are also key contributors to the chemical 
evolution of the Universe. In recent years, the relative uniformity and high luminosity of SNe 
la have made them important as "standard candles" used to make distance estimates at large 
redshifts. SNe la data are a valuable tool for measuring many cosmological parameters. The 
fact that distant SNe la appear slightly dimmer than expected for a flat, coasti ng universe 
can be used to show that the expansion rate of the universe is a ccelerating (iGarnavich 
19981 : iRiess et al.lll998al : ISchmidt. et al.lll998l : iPerlmutter et al.lll999[ ). This important result 
is combined with information from studies of the cosmic microwave background and large 
scale structure to suggest the existence of a "dark energy" in the universe. 

SNe la are not precisely standard candles, and variability in SN la observables intro- 
duces uncertainties that limit their effectiveness as distance indicators. The intrinsic disper- 
sion in the peak brightness of SNe la events can be constrained to about 0.2 magnitudes 
by calibrat ing the measured luminosity using corrections determined by the shape of the 



light curve (IPhil 



20071 : IGuv et al 



ips et al. 



2007 



1993 



Riess et al.lll995l: iHamuy et al.lll996l : iJha. Riess. &: Kirshner 



Conley et al.ll2008l ). This accuracy is sufficient for some cosmological 



analysis, but in order to achieve the level of high precision cosmology requi red for dark energy 



measurements the dispersio n must be reduced by an order of magnitude (IWeller fc Albrecht 
200ll : iKowalski et allEoOsh . 



The NIR is a productive source of information about light and intermediate mass ele- 
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ments such as He, C, O, Mg, and Mn that are undetectable or obscured by hne blending at 
other wavelengths. The progenitor of a SN la is expected to be a carbon/oxygen white dwarf 
star (C/0 WD), so that after the explosion a region where C and O are detected together 
is composed of primordial material that has not experienced nuclear burning. The extent 
of nuclear burning in the SN will consequently be defined by a boundary between a region 
in which C or O or both have been consumed and a region that contains both C and O. 
The abundance of Mn is temperature dependent so Mn can be used to probe the burning 
temperature in regions of incomplete silicon burning (IHoflich et al.lll998l . |2002| ). Si, S, Ca, 
Fe, and Co are burning products that also produce strong lines in the NIR, making data 
from this wavelength region a rich source of information about the physical characteristics 
of SNe la. 

Many of these important explosion products are revealed by spectral observations of 
SNe la made near or before the date of maximum brightness in the V-band (Vmax)- Our 
sample is an excellent resource for investigation of SNe la because it includes eleven spectra 
obtained before Vmax, and 20 spectra obtained within one week of maximum brightness. 

Supernovae expand and cool during the first weeks after the explosion which is coincident 
with the time period covered by these observations. As the envelope thins and the dominance 
of Thomson scattering in the continuum diminishes, progressively deeper regions of the 
supernova are revealed. The data show that NIR spectral features evolve consistently for 
normal SNe la, so a time sequence of spectra can be used to trace the chemical structure 
from the outer layers toward the center. NIR observations are particularly effective for this 
purpose because the optical depth is less in the NIR than at shorter wavelengths so that a 
greater radial depth can be probed with each spectrum. 

With a sample of this size, NIR data begin to realize their potential to further constrain 
estimates for the intrinsic brightness of individual SNe la by providing direct evidence of 
their physical properties. This is an essential part of the interaction between theoretical 
modeling of supernova physics and observational analysis. NIR spectra help define the 
chemical structure and constrain abundances of some burning products in SNe la. 

Relatively few NIR spectra are found in the literature and they were obtained from even 
fewer individual SNe la. S ix spectra from -8 to +8 days were obtained from SN1994D in 



narrow wavelength bands (IMeikle et a. 



19961). An excellen t time sequence for SN 1999ee 



from -9 to +42 days was obtained by iHamuy. et al.l (120021 ). but without comparison with 
other objects, the data provide no information about the typical behavior for SNe la. A fine 
set of NIR spectra w ere obtained by Gerardy from SNe la SN 1999by between -4 to +28 days 
( IHoflich et al.ll2002l ). but this object was significantly sub-luminous and does not contribute 
to the discussion of normal behavior in SNe la. Our program previously published 13 spectra 
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from 12 objects (spectra that are included in the curre nt sample) that pro vided clues for 



many of the conclusions derived from the larger sample (IMarion et al.ll2003l ). 



The spectra in this sample also provide an opportunity to search for secondary rela- 
tionships between NIR observables and other data libraries that may further constrain SN 
la physics and improve the effectiveness of SNe la as cosmological distance indicators. For 
example: this sample may provide a link to accurate "K corrections" for photometric obser- 
vations in the J, H, and K bands. 

Details of the data acquisition and reduction procedures are presented in ^ including 
the use of the Fourier transform to smooth the data. The spectra are presented in ^ 
We identify many spectral features, measure the Doppler velocities for individual ions, and 
discuss the radial distribution of explosion products. In §l]we discuss some of the physical 
implications of our results. We examine the evolution of spectral features and compare 
spectra obtained at similar ages noting the similarities and differences. Generally accepted 
characteristics of the physical properties and behaviors of SNe la are reviewed in ^to provide 
a basis for analysis of the spectra. §H] summarizes the results. 

Appendix El provides a detailed discussion of individual features in the spectra. Included 
with the appendix are Tables [5] and [6] that have estimated lines strengths of the strongest 
lines for many of the ions expected to be present in NIR spectra from SNe la. The line 
strength estimates are separately computed at two temperatures: 5,000K, which is expected 
to be reasonable for extended line-forming regions and 10,000K which is expected where lines 
are formed closer to the photosphere. 



2. Acquisition and Reduction of Data 

We obtained low and medium resolution NIR spectra from SNe la using the 3.0 meter 
telescope at t he NASA Infrared T elescope Facility (IRTF) with the SpeX medium-resolution 



spectrograph (iRayner et al.ll2003l ). The SpeX instrument provides single exposure coverage 
of the wavelength region from 0.8 — 2.5 fim. For SpeX observations using a grating and prism 
cross- dispersers (SXD mode), the average instrumental spectral resolution {Rx = A/AA) is 
750-2000 and with a single prism (LRS mode) the resolution 120-300. Most of our data 
were obtained using one of three settings: the 0.5" slit in LRS {R\ ~ 200, which makes 
the velocity resolution R^^i = (AA x c)/A ^ 750km s"^), the 0.8" slit in SXD {Rx ^ 750, 
Ryei ~ 400km s-i), and the 0.5" slit in SXD {Rx ^ 1200, R^ei ~ 250km s'^). The velocity 
resolution estimates are calculated for 1.0 fim and assume that for the LRS data we can 
accurately interpolate midway between data points (see Tables [2] and [3]) . This resolution in 
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velocity space is sufficient because the features we want to investigate have widths of several 
thousand km s~^. The LRS observing mode can extend coverage at the blue end of the NIR 
to about 0.65 fim. SpeX also contains an infrared slit- viewer /guider covering a 60x60 arcsec 
field-of-view at 0.12 arcsec/pixel. The detectors are a Raytheon 1024x1024 InSb array in the 
spectrograph and a Raytheon 512x512 InSb array in the infrared slit-viewer. 

For most observations the slit was aligned to the parallactic angle at the beginning of 
the observing sequence. In a few cases, the proximity of the SN to the host galaxy made it 
impossible to achieve the optimal parallactic alignment. However, most of our observations 
were made at less than 1.5 air-masses and the features in which we are interested are near 
the guiding wavelength so the errors due to atmospheric refraction are not significant. We 
do not attempt to compare line-strengths and continuum levels with a precision that would 
be affected by light losses due to mis-alignment of the slit. 

For bright targets, the SpeX guider effectively maintained the centroid of the target in 
the slit by guiding on the spill-over flux from the object in the slit. When the target did not 
clearly appear outside the slit, we used the SpeX guider in offset mode with another object 
in the field of the SpeX imager. Another guiding option was the IRTF optical guider and 
this was used in a few instances. 

Saturation was not a concern due to the faintness of our objects, but OH lines are 
numerous and highly variable in the NIR. To avoid an increase in background noise due to 
poor OH removal, individual exposure times are capped at 150s. A typical set is limited to ten 
exposures for a total of 25 minutes integration time in order to remain within time scales of 
atmospheric variability. The source is nodded along the slit using an A_B_B_A_A_B_B_A_A_B 
pattern. Calibration images are obtained from AOV standard stars. The standards were 
selected to be as near as possible to the time and airmass of the SN. Each observation 
set also includes calibration images from internal lamps for fiat fielding and wavelength 
calibrations. 

The data were reduced using a packag e of IDL routines specifica lly designed for the 



reduction of SpeX data (Spextool v. 3.2; ICushing. Vacca. fc Raynen 2004). These rou- 



tines perform pair subtraction, fiat-fielding, aperture definition, spectral tracing and extrac- 
tion, residual sky subtraction, host galaxy subtraction, and wavelength calibration for data 
acquired in both the prism mode and the cross-dispersed mode. Corrections for telluric 
absorption were performed using t he extracted spectrum of an AOV star and a specially 



designed IDL package developed by IVacca. Cushing fc Raynen (120031 ). These routines gen- 



erate a telluric correction spectrum by comparing the spectrum of an AOV star to a model 
AOV spectrum that has been scaled to the observed magnitude, smoothed to the observed 
resolution and shifted to the observed radial velocity. The telluric correction spectrum is 



- 6 - 



then shifted to ahgn the telluric absorption features seen in the SN spectrum and divided 
into the target spectrum. 

Spextool removes problems with background irregularities by defining the aperture 
width and background level for each subtracted image. The flux level along the slit in- 
cludes a positive image of the SN at the A position and a negative image of the SN at the 
B position. The background level in the slit is defined by Spextool using a fit through four 
regions, in two pairs, that bracket each aperture. Typically this is a linear fit, but for noisy 
spectra a second or third order fit is more effective. 

Noise levels are calculated and recorded by Spextool for each wavelength bin. Because 
the data are subject to bad pixels which can skew the combined mean value and the error 
on the mean, we substitute the median and median absolute deviation as robust estimates 
of these statistics. The median absolute deviation (MAD) is defined as: 

MAD = 1.4826 x median{\{fx - (/a)i/2)|) 

where / = flux, (/a) 1/2 is the median of the combined values, and 1.4826 is a constant 
introduced on the assumption of a Gaussian distribution of initial values. 

The error on the median at each wavelength is equal to the one-sigma noise level and 
equal to the median absolute deviation divided by Vn, where N is the number of spectra 
that have been combined. 

The data have also been cleaned by replacing spikes greater than twice the local noise 
level (Figured]). The replacement value from removed spikes is the mean of 40 adjacent data 
points for SXD and 8 adjacent points for LRS spectra. 

Observing time scheduled months in advance presents unusual problems when working 
with transient objects such as SNe. We recognized that the number of potential targets would 
be improved by more frequent but shorter observing runs. Working with IRTF personnel we 
employed a remote observing program to obtain data on systematically scheduled observing 
dates at a frequency of 8-12 times per year for five years. We obtained high quality NIR 
spectra from whatever SNe la were available at the scheduled time. In most cases, each 
target was accessible only once during our observing dates. As a result, we obtained a single 
epoch or "snapshot' from most of our targets. 
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2.1. Using the Fourier Transform technique to smooth spectra 

Using a Fourier transform (FT) permits us to smooth the data by removing the highest 
frequency portion of the signal, leaving only the lower frequency components which are 
the actual information from the target. Since the features we are studying are very broad 
compared to the spectral resolution, this is an excellent way of reducing the noise to allow 
more consistent measurements and more effective comparisons between spectra. 

The energy flux from SNe la drops off rapidly with increasing wavelength, and as a 
result the signal-to-noise in our data diminishes at longer wavelengths. Due to this different 
quality of the data, we use different parameters for the FT smoothing in different parts of 
the spectra. SXD spectra contain ^ 3,900 data points covering the wavelength range 0.8-2.5 
Hm. This large number of bins allows us to separate SXD spectra into three regions: 0.7-1.35 
/xm, 1.35-1.8 yum, and 1.8-2.5 /im. LRS spectra span the same wavelength band with ^ 540 
data points. We separate LRS spectra into two sections with the division at about 1.35 /im. 

Each region of a spectrum is individually flt with a straight line between the endpoints 
and normalized to a flat continuum with the continuum equal to zero. Setting both ends of 
the spectrum to zero eliminates most of the ringing that can be introduced by the transform 
procedure. After normalization, a fast Fourier transform routine is used to change the 
spectrum to frequency space. 

The power spectrum (signal power by frequency) is obtained by taking the square of 
the absolute value of the frequency spectrum. A log plot of the power spectrum reveals the 
noise level forming a horizontal line through the central region of the spectrum. The signal 
power as a function of frequency descends rapidly from a maximum power in the lowest 
frequencies located at the edge of the spectrum and it crosses the noise level usually within 
the flrst 50 frequency bins. Figure [2^ shows a log plot of the 50 data points (out of 851) 
in the power spectrum derived from the region 0.8-1.35 /im in the SXD spectrum obtained 
from SN 2005am at -4 days. The noise level (N) is taken to be constant with a value equal to 
the mean of the 100 data points immediately to the right of the region displayed (50-150). 
N is marked on the flgure by a horizontal dashed line. 

We estimate that the slope of the signal for this data intersects the noise level near 
frequency bin 27. We represent the signal level (S) for each frequency bin using a straight 
line in log space from the y-intercept (~ —4.0) through the intersection point: (27, -12.5). To 
avoid an abrupt transition from the signal line to the noise at the intersection, we generate 
a filter using S/(S-|-N) and multiply the frequency spectrum at each bin by the filter at that 
bin (Figure[2]D). We explored using a Wiener filter (S"^/ (A^^ -|- S"^)) and found no measurable 
difference in the resulting smoothed spectrum. The result of applying the filter to the FT 
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data is shown in Figure [3] with the original frequency spectrum plotted in black and the 
filtered spectrum is superimposed in red. 

As a test of the consistency of our method, we compare the effects of different smoothing 
parameters on differently shaped absorption features from the same data (SN 2005am at - 
4d). The smoothing parameter is defined as the number of frequency bins to be included in 
the signal (S). Figure H] shows raw data plotted in black with the results of various smoothing 
parameters superimposed in colors. The top panel displays the feature from Mg II at 0.9227 
fim which is asymmetric due to impingement on the red side from emission due to the Ca II 
IR-triplet. The bottom panel shows the feature from the Mg II 1.0927 fim line which is more 
nearly symmetrical. 

The FT smoothed spectrum, derived according to our standard procedure described 
above and using 27 bins, is plotted in black. We reduce the number of bins until the 
smoothing obviously compromises the positions of the absorption minima. Different colors 
are used in the figure to represent the spectrum smoothed using different parameters. Table 
[T] gives the measured velocities with the different smoothing parameters. The table shows 
that a very large change in the smoothing parameter, beyond any reasonable choice made 
by the user, is required to change the measured velocity by more than a 200 km s~^. 

When the smoothing parameter is increased, the change in measured velocity is very 
small. The primary effect of too large a value of the parameter is that local noise may be 
included in the smoothed spectrum. However the features we are measuring are much larger 
than any local spikes. Consequently, when there were ambiguities in the choice of smoothing 
parameters, we were careful to avoid choosing smaller values. Interpretations made by 
individual users may vary the total number of bins selected by 3-5 in the smaller direction 
and 5-10 in the greater direction but the figure and table clearly show that differences on 
that scale have no significant effect on the results. An exaggerated change on the order of 
40% (from 27 to 16 bins) is far beyond the uncertainty in our parameter choices and it only 
moves the measured velocity by 100 km s~^ for the more symmetrical feature in the bottom 
of Figure m and 900 km s^^ for the less symmetrical feature in the top of Figure HI 

The amount that the spectral resolution is degraded by smoothing depends on the 
original resolution and the amount of smoothing required. In the example here with the 
data from SN 2005am at -4d and a smoothing parameter of 27, the instrumental resolution 
of ~ 250 km s~^ is reduced to ~ 1500 km s~^, but the measured location of the absorption 
minimum only changes by 200 — 300 km s~^. 

The reverse Fourier transform of the filtered frequency spectrum returns a smoothed 
data spectrum that greatly improves our ability to interpret the features. The top spectrum 



- 9 - 



in Figure [5] shows the raw data from a relatively high quality spectrum obtained from SN 
2005am at -4d plotted in green, with the smoothed spectrum plotted in black over the 
original. The lower spectrum in the figure shows data from a noisier spectrum, SN 2004bl 
at -2d, and again the smoothed spectrum is plotted in black over the original data. The 
figure illustrates the effectiveness of this technique to facilitate inspection and comparison 
of spectra in the sample. 

2.2. Fitting a local continuum 

The energy flux from SNe la decreases rapidly with increasing wavelength. Within 
this steeply sloping continuum, the wavelength of minimum flux may not correspond to the 
wavelength of the maximum departure from the continuum. For absorption features, this 
effect is not significant where the continuum is relatively fiat, as in the optical region of the 
Si II line at 6355 /im, but it can potentially affect features in the NIR by ~ 0.005 /im or 
~ 1,500 km s^^. To remove this influence, we estimate the location of the continuum in 
the region of each absorption feature and normalize the spectrum to a flat continuum before 
identifying the wavelength of absorption minimum. 

Line-forming regions in SNe la are extended in velocity space ( §5.3p . Thus a full analysis 
requires measurement of the full range of Doppler velocities between the absorption minimum 
and the detection limit of the blue wings for each feature. Due to uncertainties introduced 
by the sloping continuum, this region can be more accurately measured by normalizing the 
data to a level continuum. 

Figure illustrates the difficulty of fitting a continuum to SNe la spectra over the wave- 
length region 0.8-2.5 fim. Four spectra from different epochs are displayed with estimated 
continua superimposed in red (-3.0 power law) and green (-4.5 power law). It is clear that a 
single power law does not fit any spectrum over the entire region. Although we have been 
as consistent as possible with our observing and reduction techniques, it is possible that 
some of the differences in slope and shape between spectra may have been introduced by 
observation or data processing choices. There are however, obvious changes in the overall 
shape of the spectra with time and there appear to be measurable differences in the shapes 
of spectra from different objects at similar epochs. 

For the purpose of comparing features from the same line in different spectra, we aban- 
don our efforts to fit a continuum over a large wavelength range and concentrate on fitting 
a local continuum in the vicinity of each feature. Figure [7] shows the same spectra as the 
top two panels in Figure E] with local continua fit to the spectra in the wavelength regions 
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close to each absorption feature. This technique produces better results than fitting a single 
continuum to the entire spectrum. The fitted regions are expanded in Figure [8] 

There remains ambiguity about how to fit the local continua because of differences 
in the spectra within each region but variations within the limits of our estimates do not 
significantly affect the results. For example: we assume that the peaks near 0.85 and 0.93 fim 
rise above the continuum as emission components of P-Cygni profiles from the very strong 
Ca II triplet (mean=0.8579 ^m) and the strong Mg II line (0.9227 /xm); however, it is not 
obvious where the true continuum level lies. Our estimate for the location of the continuum 
in this region appears as the green line in both panels of Figure [71 

The absorption features in Figure [7]with minima near 1.05 and 1.20 fim display different 
profiles in the two spectra. It is not clear how much of the bump on the red side of the 
absorption is emission and how much is continuum. The estimated locations of the local 
continua for these regions are shown as red and blue lines in Figure [71 

Figure [8] expands the spectral regions with locally fit continua from Figure [71 but with 
the continuum normalized to zero. All figures on the left are from the SXD spectrum of SN 
2002en at -4 days. The figures on the right are from a LRS spectrum of SN 2002cr at -5d 
days. The top row corresponds to the regions marked in green on Figure [71 the middle row to 
the red regions and the bottom row to the blue. The FT smoothed spectrum is superimposed 
in red on the spectra in all figures. We tested this method by gradually changing the fitted 
continua until it was clearly beyond a reasonable fit to the data. The largest deviations in 
the wavelengths of the measured minima due to the manual placement of the continua was 
found to be ±300 km s~^ and for most features the deviations are smaller. 



3. The Spectra 

Our sample consists of forty-one (41) NIR spectra from SNe la, obtained between four- 
teen days before (-14d) and seventy-five days (+75d) after Vmax {FigmeMj- All spectra in our 
sample were obtained at the IRTF using the SpeX instrument. Observation and reduction 
practices are discussed in ^ 

The spectra are listed in Table [H by epoch with respect to Vmax and Table [31 in order 
of discovery. The Tables include: SN names, estimated ages at observation (with respect 
to Vmax), the dates of observation, the average instrumental resolution observing mode, the 
estimated signal-to-noise ratio for wavelength regions approximately corresponding to Y, J, 
H, and K bands, and the redshift (z) of the host galaxy. The epoch of observation with 
respect to Vmax is estimated by a variety of means and we have attempted to provide as 
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accurate information as possible. More complete light curve data for many of the SNe in 
our sample are in preparation and the observation date for each spectrum is included in the 
tables in order to facilitate comparison when these data become available. 

In all figures containing multiple spectra, we arrange the spectra in sequence with the 
earliest spectrum at the top according to the time of observation with respect to Vmax- The 
abscissa is wavelength, in microns, shifted to the rest frame of the host galaxy and the 
ordinate is log flux. 

Each spectrum has been normalized to 1.0 at 1.0 /im. The names of the SNe and their 
age when observed (in days relative to Vmax) are displayed to the right of each spectrum in 
the same color as the spectrum. The colors for each spectrum are consistent throughout the 
colored figures to facilitate identification. The black line superimposed on the original data 
is the smoothed spectrum obtained by Fourier transform as described in §21 

Atmospheric opacity depends strongly on wavelength in the NIR. A transmission spec- 
trum is provided at the top of some figures to help identify regions where opacity reduces 
the signal. Complete transparency is at the top axis of the plot and full opacity at one unit 
down. The transmission spectrum has been shifted using z = 0.02 which is the approximate 
mean redshift of the SNe la in our sample. The noise levels in the data are very high through 
regions where the opacity exceeds 20% (~ 1.30 — 1.42 /im and 1.78 — 1.90 //m). We have 
removed the data in these regions from all spectra. 

Order transitions in SXD spectra occur near 0.95, 1.12, and 1.45 fim. Small local 
discontinuities are possible at these locations where the orders are merged. This problem is 
more pronounced in the noisier spectra. SXD spectra also have gaps between ^ 1.35 — 1.42 
where there are no data between orders. The gap locations were selected to correspond to 
regions of high atmospheric opacity. 

Type la spectra are characterized by the steadily diminishing flux at longer wavelengths, 
so the S/N ratio degrades with increasing wavelength as the level of the signal approaches 
background noise levels (Tables [3] and [2]). 

3.1. Line Identification by Epoch 

For spectral analysis we removed seven spectra that have S/N ratio less than 10 in 
the region 1.10 — 1.30 fim from the complete sample (see Table [2]). This wavelength region 
approximately corresponds to the J-band and contains many of the most important features 
under discussion. These seven spectra are displayed using the color yellow in Figure M The 
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thirty-four remaining spectra are presented in Figure [TOl 



For analysis and line identification, we divide the remaining spectra into three groups 
selected by epoch and spectral features. The boundaries between these phases is not always 
distinct and there is some temporal overlap between the phases. We define the groups as: 
Photospheric Phase (Outer Layers of the SN; 18 spectra obtained between -14d and +5d), 
Extended Photospheric Phase (Intermediate Layers of the SN; 10 spectra between +4d and 
+24d), and Transitional Phase (Inner Region of the SN; 7 spectra obtained between -|-37d 
and +75d). 

We identify many of the prominent features, measure their Doppler velocities, and 
provide a synopsis of spectral development in each group. Appendix |X] provides a more 
detailed analysis of the spectra and features with possible alternative line identifications and 
a discussion of each ion present and the physical implications. We display the spectra in log 
space to assist with the identification of features at longer wavelengths, but care must be 
taken when comparing the relative strengths of features from widely separated wavelength 
regions. 

Figure [TT] shows the spectrum from SN 2002fk at -14d, plotted in linear space and also 
in log space. Spectral features near 1.6 and 2.1 /xm are clearly smaller than the features 
appearing at 1.05 and 1.22 /im when plotted in linear space, but these features appear to be 
larger when plotted in log space. 

Acceleration of matter in the explosion lasts for only a few seconds but withi n minutes 
the envelope expands by several orders o f magnitude from the initial WD radius (IKhokhlov 
199ll : iHoflich. Mueller fc Khokhlovl Il993l ). The geometry of the expanding SN is assumed 
to be approximately spherical because the progenitor WD was nearly spherical, and the 
ejecta continue to move in hom ologous expansion with radial velocities proportional t o the 
radial distances from the center (IKhokhlov. Oran fc Wheeled 1 199 71 : iGamezo et al.ll2003l ) and 
see discussion in ||5l A continuous range of expansion velocities can therefore be used to 
represent contiguous physical space in the supernova, and we often refer to comparative 
radial distances by their measured velocities. For example, a layer, or shell, of material in 
the spherical envelope may be described as occupying a region in space between 10,000 and 
12,000 km s-\ 



Polarization studies indicate that sorne SNe la, at least in th e outer layers, may not 
be completely spherical (IWang. et al.ll2003l : IWang fc Wheeled 120081 ). Possible explanations 
for the asphericity include rotation effects of the progenitor, mergers of two C/0 WD stars, 
high-velocity filaments of explosion products, and influences on the ejecta by the binary 
companion. The observed polarization appears to be at maximum in high-velocity Ca II 
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regions. The fact that high-velocity Ca II components are present in most of the early spectra 
in our sample suggests that NIR data will be valuable for further study of asymmetry in 
SNe la. 

After a preliminary identification of an association between a spectral feature and an 
individual atomic line, it is prudent to examine the spectrum for other lines from the same 
ion. Comparison of relative line strengths can establish whether other lines from that ion 
would be expected to produce detectable features and confirm the initial identification. To 
easily make these comparisons we have calculated line strength estimates for several of the 
strongest lines at two temperatures: 5,000K, which represents the temperatures found in 
extended line- forming regions, and also at 10,000K which represents temperatures found 
closer to the photosphere. These temperatures are sufficiently separated to explore the 
effect of temperature on the relative strengths of lines from the same ion. It is not our 
intention to precisely define temperatures in the line-forming regions. 

The relative line strength estimate uses the LTE approximation for the fraction of ions 
from a population that will be excited to the n*'* level. The population of level n is given by 
the Boltzmann factor: A^„ oc * lO"^"/'^"^, where /„ is the oscillator strength, is the 
statistical weight, Xn is the excitation potential, and T is the temperature. 

In Table |5] the estimated line strengths for ions is ordered by strength at 5,000K and 
in Table El they are ordered by strength at 10,000K. Computed values for different lines 
from each ion are normalized to the strongest line and rest wavelengths are given in air. In 
situations where multiple lines from the same ion are close in wavelength, we calculate a mean 
value, weighted by oscillator strength, for two or more lines. For these blends, the \og{gf) 
value and excitation minima are listed for the strongest line in the blend. Our generalized 
treatment of line strength estimates has limitations but it is not model dependent. The 
use of NLTE departure coefficients introduces uncertainty in the estimated line strengths by 
factors of 2-10. Uncertainties of this magnitude do not change our results. 

Transitions from one ionization stage to another are expected to be abrupt in both time 
and velocity space due to the dominance of radiative transitions over collisional processes. 
Fe-group lines dominate the photo-ionization spectrum so ionization boundaries for other 
elements are influenced by flux levels from the Fe transitions. Co nsequently most li ne- 



forming regions will contain only one ionization stage for each atom (IMarion et al.ll2006l ). 



-14- 



3.1.1. A Very Early Spectrum (-14d) 

Figure [11] shows the spectrum of SN 2002fk, obtained 14 days before Vmax which is only 
five or six days after the explosion. This is a single spectrum obtained in SXD mode. For a 
detailed discussion of individual features, please reference §A.1[ 

Probable Identifications: Mg II, Si II, and Ca II. 

Possible Identifications: O I, Si III, Fe II, Fe III, and Co II 

Synopsis: The measured Doppler velocities in this spectrum are 10,000-12,000 km s~^ 
which is slower than the measured velocities for spectra from other SNe la obtained at 
later epochs. There is also possible evidence for multiple ionization stages of silicon. These 
unusual characteristics may be caused by lower temperatures within the envelope at this early 
epoch than will be found a few days later. The lower temperatures allow the photosphere 
to temporarily recede more deeply into the ejecta. Soon after the epoch of this spectrum, 
energy deposition from radioactive decay of ^^Ni increases the temperature and opacity which 
pushes the photosphere out to higher velocities (A more detailed discussion of this effect is 
found in §5.31) . 

He I and C I are not detected in this spectrum. O I is possible from the 0.9264 fim line 
at a Doppler velocity of 11,500 km s~^ in a blend with Mg II. Mg II is clearly detected from 
lines at 0.9227 and 1.0927 /im and is probable from lines at 1.0092 and 1.6787 /im. Doppler 
velocities for Mg II are 10,300-11,600 km s~^. A feature due to Si II from the 0.9413 fim line 
is likely to be present at 11,300 km s~^. Si III is possible but at a low Doppler velocity of 
8,300 km s~^. The Ca II IR-triplet is strong at 12,300 km s~^. Fe II is unlikely and Fe III is 
possible but at a low Doppler velocity of 8,100 km s~^. Co II is detected at 12,000 km s~^. 

3.1.2. The Photospheric Phase (Outer Layers; -l^d to +5d) 

Figure [T^ shows eighteen spectra obtained between -14 and +5 days relative to Vmax- 
During this epoch the photosphere is expected to be in the outer layers of the SN and the 
line-forming regions will be close to the photosphere. Mg II is prominent from multiple lines 
in all spectra in this group. Six of the spectra were obtained in LRS mode which extends 
blueward to about 0.7 /im at R=200 (compared with R=1200 for the SXD mode). Figure 
[T^ shows four representative spectra from this group: SNe 2004bw (LRS obtained at -9d), 
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2005am (SXD obtained at -4d), 2000dm (SXD obtained at +0d), and 2001bf (SXD obtained 
at +4d). Only the region 0.70 — 1.32 /im is displayed in Figures [12] and [13] in order to enhance 
details in the area where most of the features are under discussion. Features are labeled in 
the figure for which we have reasonably confident identifications. For a detailed discussion 
of individual features from Photospheric Phase group, please reference §A.2[ 

Velocity references are compiled in Table H] 
Probable Identifications: O I, Mg II, Si II, Ca II 
Possible Identifications: Mn II, Fe II, Fe III, and Co II. 

Synopsis: Doppler velocities for lines identified in this group are between 10,000-14,000 
km s~^ with lower velocities from both spectra of SN 2002cr (-7d and -5d) and from SN 
2001bf (+4d). 

C I is not detected in any of the spectra ( §4.3p . Strong features near 0.75 /xm due to 
O I are detected in all spectra in this group that include the wavelength region. Mg II 
produces multiple features in all spectra from this group and the Mg II lines clearly correlate 
by velocity and strength. The Ca II IR triplet is prominent in all spectra from this group 
and a high-velocity component of the Ca II feature is present to some extent in all spectra 
that include the wavelength region shorter than 0.8 /im. Si II is detected either alone or in 
a blend with Mn II. Fe II is likely in blends with other atoms, as is Co II. 

A bump appearing near 1.25 fim is strongly peaked in some of the spectra, unlike emis- 
sion from a P-Cygni profile. This feature does not have a confident identification although 
it is observed in all spectra throughout our entire sample from -14d to +75d. The feature 
near 0.97 /im becomes stronger after maximum light due to contributions from Fe II. An- 
other feature of uncertain origin appears at the same time near 1.48 /im. Beginning about 
maximum brightness, the spectra begin to be affected by changes in the apparent continuum 
due to line blanketing from Fe-group lines. "Emission" bumps begin to appear near 1.25, 
1.52, and from 1.6 - 2.1 /zm and in broader regions (see discussion §5.3p . 

3.1.3. The Extended Photospheric Phase (Intermediate Layers; +4d to +24d) 

Figure [H] shows ten spectra from SNe la obtained between 4 and 24 days after Vmax- 
During the phase represented by this group, the shape of the continuum becomes influenced 
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by line-blanketing from Fe-group lines which extends the effective photospheric radius (see 
§5.3p . Spectra obtained in LRS mode extend blueward of 0.80 fim and all others were 
obtained in SXD mode. Figure [T^ shows spectra from SNe 2000dk (+5d), 2001bg (+10d), 
SN 2004ab (+18d), and SN 2004da (+24d) as representative of this group with features 
labeled for which we have confident identifications. For a detailed discussion of individual 
features from this group, please reference §A.3[ 

Velocity references are compiled in Table |H 
Probable Identifications: O I, Mg II, Si II, S I, Ca II, Fe II, and Co II. 
Possible Identifications: Si I, Mn II 

Synopsis: During the temporal period represented by this group, transitions are made 
from spectra dominated by O I and Mg II lines to spectra with lines from S I, Fe II, and 
possibly Mn II and Si I. Measured Doppler velocities remain above 9,000 km s~^ for most 
ions with the exception of Fe II which begins near 9,000 km s~^ at earlier epochs but slows 
to less than 3,000 km s"^ by +18d. 

O I continues to be detected up to +10d and is possibly identified to +18d. Mg II is 
likely up to +5d but is absent thereafter. The absorption near 0.89 fim shifts to 0.90 and 
apparently makes a transition from Mg II to S I with possible contributions from O I. This 
feature becomes much weaker after +5d which is consistent with losing the signal from Mg II 
and O I. 

Si II and Mn II are possible alone or in a blend producing the feature found near 0.92 
fim that is not present in all spectra. The Ca II IR triplet remains prominent and the data 
continue to separately resolve the two strongest components of this feature. Fe II is detected 
from the line at 0.9998 /im in all spectra through +18d, and Co II is strong in the spectra 
from this group obtained at +14d and later. 

Individual Fe-group lines are complicated by many potential transitions at the same 
wavelength and nearly the same strength within each multiplet so the simple formula we 
have been using to estimate line strength may not be accurate for Fe II as it was for smaller 
atoms. In the spectra from this group, the 0.9998 fim line is the only Fe II line detected (+4d 
through -|-18d), although the calculated line strengths suggest that lines at 0.7462, 0.7712, 
and 1.0500 fim should also be detected ( §A.3p . 

During the Extended Photospheric Phase, the shape of the continuum becomes influ- 
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enced by line-blanketing from Fe-group lines (see §5.31) . Looking at this group in the context 
of the entire sample (Figure [TOj) . we can see how unresolved groups of Fe lines change the 
effective shape of the continuum during the epoch covered by these spectra. Significant peaks 
arise near 1.25, 1.55, 1.75 2.0, 2.1, and 2.25 fim. Thus many of the apparent "absorption" 
features may be interpreted as gaps between regions of inflating continuum. For example: 
minima found near 1.12, 1.21, 1.50, and 1.65 /xm are probably formed by local flux increasing 
on both sides of an apparent minimum rather than by absorption at that location. 

Another example can be seen in Figure fTOl where the entire region from 0.9 to 1.3 /xm 
displays a distinct change from a concave bend below a straight line in the epoch +0d to 
+4d, to a nearly straight pattern from +4d to +10d, and finally to a convex bulge above a 
straight line after +10d that continues in later epochs. 

This group exhibits a reduced continuum level from 1.2 — 1.5 fim as compared with the 
earher spectra. Figures [H [TOl [HI and [15] clearly show this effect. This depressed continuum 
level is due to the region having fewer blends of iron group lines compared to adjacent 



wavelengths (j Wheeler et al .1119981 ). Note that the spectrum of SN 2001bf (+4d) does not 
exhibit any deficit in this region, even though the epoch is estimated to be only one day 
from that of SN 2000dk (+5d) in which the effect is pronounced. 

We note that the notch found at 1.65 fim between two emission peaks, occurs at ap- 
proximately the same wavelength as the Mg/Si/Co absorption in earlier spectra, but the 
features are probably not related. The spectrum of SN 2000dk at +5d is the oldest to show 
the strong notch indicative of Mg II absorption at this location and later spectra have a 
smoother, shallower depression in this wavelength region. The spectrum from SN 2000dk at 
+5d is also the oldest spectrum to show Mg II absorption from lines at 0.9227 and 1.0927 
yum (see Figure [T^ . 

The steep slope on the blue side of the peak found near 1.55 fim defines the transition 
from partial to complete silicon burning. Due to the high optical depth, this feature is formed 
close to the outer edge of the Fe/Co/Ni core at a rest wavelength of ~ 1.57 fim. Precise 
determination of the transition feature is not possible but velocities measured to a point 
approximately half way up the slope are found between 9,900-12,100 km s~^ for spectra in 
this group obtained from -t-4 to -t-18 days after maximum light. 



3.1.4- Transitional Phase (Inner Region; +37d to +75d) 



Figure [16] shows seven spectra from our sample obtained between 37 and 75 days after 
Vmax- This group displays a lot of new structure relative to earlier NIR spectra and the 
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features are remarkably similar from one spectrum to the next. The spectrum of SN 2004bv 
at +60d fits well with the others in this group although the spectrum from this SN at -5d 
was unusual (see Figure [T^. Spectra obtained in LRS mode extend blueward of 0.80 /xm. 
All others were obtained in SXD mode. For a detailed discussion of individual features from 
this group, please reference §A.4[ 



Synopsis This epoch is a transitional phase between the photospheric era and a true 
nebular phase. This is the same time frame in which unreddened SNe la have a very low 
dispersion for B-V colors (Lira, 1995). The common assumption about this era is that 
expansion has reduced the opacity beyond the point where a distinct photosphere can be 
defined and individ ual absorption features can no longer be produced. Recent work by 



Branch, et al.l (120081 ) suggests spectra from this epoch can be modeled by absorption features 



from permitted lines. 

The apparent emission peaks at 1.60, 1.75, 2.15, 2.25, and 2.35 /im are the result of 
increased flux due to line blending from iron group elements as explained in §5.31 The 
radius of optical depth unity has increased at these wavelengths producing a larger effective 
radiation area and an increased flux. 

Fe II lines at 0.7712, 0.9998, and 1.0863 fim appear in absorption beginning about 30 
days after maximum light. In spectra from our sample obtained before +37d, the 0.9998 /im 
line is the only Fe II line detected, although the calculated line strengths suggest that lines 
at 0.7462, 0.7712, 1.0500, and 1.0863 fim lines should also be detected. 

Co II lines are clearly detected from lines at 1.7772, 1.7462, 2.2205, 2.4596, and 2.3613 
/xm at velocities near 9,000 km s"^. Co II is found in wavelength regions where a pseudo 
photosphere is created at extended radii (1.5-1.85 and 2.05-2.45 /xm). Strong absorptions 
are detected from several other Co II lines. 

Figure [T71 shows evidence for Co II and Fe II lines in 13 spectra from our sample obtained 
between +14d and +75d with respect to Vmax- The figure includes two a dditional spectra 



from SNe 1998aq and 1998bu (plotted in yellow) that were obtained by iRudy fc Puetter 
(private communication) and have not been previously published. These spectra are added 
to verify the consistency of our results. The positions of strong Co II and Fe II lines drawn 
on the figure to indicate 11,000 km s~^ at the top spectrum and 6,000 km s~^ for the bottom 
spectrum. This corresponds to the reduction in measured velocities during this epoch. 

The steep increase in flux found near 1.55 /im is used by some models to define the 
transition from partial to complete silicon burning. This feature is well defined up to about 
+60d. As shown in Figure [161 the spectrum from SN 2004bv at +60 displays this feature 
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but it is smeared out in the spectrum of SN 2001fe, also obtained at +60d. 



4. Results 

4.1. Uniformity of Spectral Evolution 

Figures [9], [THl [12l [HI and [16] reveal that the development of spectral features with time 
is very consistent throughout the sample. The data show that spectra obtained from normal 
SNe la at close to the same age have similar but not identical characteristics. There are 
many small variations in features between adjacent spectra and a few instances where one 
spectrum displays a prominent feature that is absent or weak in other spectra with similar 
relative ages. Nevertheless, within uncertainties introduced by variations in data quality 
and dating methods, the progress of feature development and overall shape is remarkably 
coherent from spectrum to spectrum in our sample. That means that NIR spectra from 
different objects can be directly compared and arranged sequentially to mimic a time series. 
Even spectroscopically unusual SNe la such as 2004bv and 2004da, which are missing some 
of the features found in most other SNe la, fit well with the general profile of other spectra 
obtained at similar epochs. The agreement in spectral features from so many different SNe 
la eliminates the possibility that observations of finite line-forming regions ( §4.21) are due to 
clumps of material. 

We test this hypothesis by comparing our "sna pshot" spectra f rom s eparate SNe la to 



an excellent time series of NIR spectra obtained by iHamuy. et al.l (|2002| ) from SN 1999ee. 
Five spectra from this series obtained between -lid and -|-13d are displayed with twenty-four 
spectra from our sample in Figure [181 Our data include a series of four spectra from SN 
2005am (-4d to -|-14d). The progression of feature development in a time series of spectra 
obtained from individual SNe la matches well with the sequence of spectra obtained from 
diverse objects. It is clear that spectral development of normal SNe la in the NIR is very 
consistent. 

Although the overall NIR spectral evolution for SNe la is uniform in basic shape and 
features, there is obviously a great deal of diversity in the class of "normal" SNe la. For ex- 
ample: our most confident line identification (apart from the Ca II triplet) is Mg II at 1.0927 
fim. Yet spectra obtained within a day or two of each other sometimes exhibit dramatically 
different depths and widths for the absorption features produced by this hne (see Figures 
[12] and [13]) . Other examples of diversity include the anomalous velocities of SNe 2002cr and 
2001bf, the weak features in the -4d spectrum from SN 2004bv, the variations in size and 
shape of the notch at 1.6 fim over a range of epochs, and the vagaries of the unidentified 
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bump at 1.25 jjm. The presence of such variabihty provides motivation for further research 
and data collection efforts. Despite these variations in individual features, the overall spec- 
tral shape and the general behavior of most features permits placing a spectrum in temporal 
order within a few days relative to Vmax- 



4.2. Measured Doppler Velocities 

The Doppler shift from the rest wavelength of a line measures the velocity of the ab- 
sorbing material along the line of sight to the observer. As discussed in §5.31 when the 
photosphere is close to the line-forming region, the observed Doppler-shift of the absorp- 
tion minimum of the spectral feature corresponds to the velocity of the photosphere. At 
later times expansion moves the line-forming region away from the photosphere. If there is 
a distinct discontinuity in the composition of the line-forming region, then a boundary of 
the absorbing layer for a particular ion is clearly defined and the velocity of the absorption 
minimum for an ion is determined by the inner edge of this layer. The absorption feature 
will decrease in strength as the column density drops, but the velocity of the absorption 
minimum will no longer change until the ion is undetectable. 

Table H] and Figure [12] show the measured Doppler velocities for O I, Mg II, Si II, S I, 
Ca II, and Fe II. Distinct minimum velocities, as measured from the absorption minima 
of spectral features, exist for all of these intermediate mass ions. These velocity minima 
correspond to the location of the inner edge of a layer of line-forming matter in velocity 
space and physical space ( §5.3p . Figure [19] shows that the line-forming regions from two O I 
lines and three Mg II lines occupy similar velocity spaces at ~ 10,000-15,000 km s~^. The 
velocity spaces for Si II and S I are slightly lower but overlap the O I and Mg II region at 
~ 8,000-12,000 km s~^. The Ca II velocities are measured for the IR-triplet and may not 
be produced by fresh burning products ( §5.31) . This may account for the fact that minimum 
velocities for Ca II are not as consistent as they are for some other species. Fe II appears to 
have no minimum velocity as the detected velocity descends from about 12,000 km s~^ at 
10 days before maximum brightness to about 2,000 km s~^ at 20 days after maximum. 

The distribution of burning products and radially-stratified chemical structure is con- 
sistent with the explosion occurring in a C/0 WD with a radial density gradient and lower 
densities on the outside ( §5.2p . The data appear to eliminate the possibility of large scale 
mixing of lower mass elements to the center or Fe-group elements to the extreme outer layers. 
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4.3. Non-detection of Carbon 



Figure [20] displays tlie eigliteen spectra from Figure [12] with the locations of strong C I 
lines marked at Doppler velocities of 11,000 km s~^. The strongest lines at 5,000 K, as listed 
in Table [5] are marked in red on the figure and the weaker lines appear in black. Small 
features appear in some spectra near the expected location for lines from 0.9406, 1.1330, and 
1.1754 /im, but if these lines produce absorption features then the strongest lines at 1.0693 
and 0.9093 /im should be detected and they are not. The possible feature for the 0.9406 /im 
line is more likely attributed to Si II at 0.9413 (A complete discussion of this group of 
spectra is found in §3.1.21 and IA.2I ) 

Carbon is not detected in these data, a result that agrees with previous results wit h 



NIR data establishing very low limits on carbon abundance in SNe la (iMarion et al. 



Similar results have been obtained using C II lines in optical spectra (ITanaka. et al 



2006). 



20081 ) 



but the much higher excitation potentials for C II lines in the NIR make it unlikely that 
C II will be detectable in these NIR data. It is clear however that NIR data show that 
carbon burning products, O and Mg, are plentiful in the outer layers and that C is absent. 
The line strengths for C and O will be similar for the same ionization stages, and if O I 
in the da ta indicates unburn ed material, then C I should also be present in detectable 
amounts (IMarion et al.l 120061 ). This result suggests that the entire progenitor is burned in 
the explosion. 



5. SNe la Physics 
5.1. Generic Model for SNe la 



Many years of optical observations support the widely accepted model that SNe la are 
the result of a thermo nuclear explosion in a carbon-oxygen white dwarf star (C/0 WD) 
( iHoyle fc Fowled ll960l ). The WD mass at the time of explosion is expected to be near the 
Chandrasekhar limit {Mch ~ I.4M0). The progenitor WD begins with a mass well below 
Mch and may grow by accreting matter from a red giant or main sequence binary companion 
by Roche lobe overfiow. Alternative models for SNe la create progenitors with a He star 
or another WD as the donor companion. These evolutionary tracks will also produce a 
C/0 WD near Mch- Burning of accreted H or He on the surface of the WD conve rts it to 
C and O in approximately equal abundance (iDommguez. Hofiich &: Stranieroll200ll ). Since 
the WD is supported against gravitational collapse by the degenerate electron gas, the final 
structure does not depend on the history of the progenitor. Given the similar quantity and 
composition of materials, it is reasonable to assume that all progenitors will explode with 
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similar properties 



The primary fuel in the combustion of SNe la is carbon and oxygen which produces an 
ash that consists predominately of elements between silicon and nickel. The composition of 
the final burning products does not affect the energy production because very little additional 
binding energy is released in the final transition from Si-group to Fe-group elements. If a 
C/0 WD progenitor with mass ~ IAMq burns its entire mass to at least Si-group elements, 
^ 2 X 10^^ ergs of nuclear energy will be produced. That is sufficient energy to overcome the 
gravitational binding energy of 5-7 xlO^° ergs and provide the kinetic energy to accelerate 



the ejecta to observ ed expansion velocities (IHofiich. Mueller fc Khokhlovlll993l : lHofiiclj|l995 
Hoflich et al1l2002h . 



The results described in §4.31 suggest that the progenitor is almost entirely consumed in 
the explosion. The nearly complete burning of progenitors with very similar characteristics 
can explain why the total production of nuclear energy is approximately constant in SNe la. 



5.2. Chemical Structure of SNe la 



Optical spectra of SNe la reveal a range of elements produced by nuclear burning of 
carbon and oxygen at densities of 10^—10^ g cm~^. To agree with observations of intermediate 
mass elements, burning conditions must exist in SNe la at densities less than 10'' g cm~^. 
Since the initial density of a C/0 WD near Mch is ~ 10^ g cm~^, it is necessary for early 
burning to lift the WD in its gravitational potential causing it to expand and reduce the mean 
density ( Hofiich et al. 19981 ). The quantities and distribution of different burning products 
are determined by the relation between the hydrodynamical and nuclear time scales. The 
hydrodynamical time scale is ~ 1 second, but the nuclear rates vary with density. The range 
of densities implied by observations inclu des three distinct regions that produce different 



products from ther monuclear combustion flHoflich. Mueller fc Khokhlovlll993l : lHofiichlll995 
Hoflich et al1l2002h . 



At densities greater than 10^ g cm~^, temperatures exceed 5 x lO^K. Thermonuclear 
burning under these conditions proceeds uniformly to nuclear statistical equilibrium (NSE) 
and produces Fe-peak elements (Fe/Co/Ni). In the case that Phum > 10^ g cm~^, electron 
capture results and stable isotopes are produced. For 10^ g cm~^ > pbum > 10^ g cm~^ the 
final burning product is predominantly ^^Ni. The existence of low and intermediate mass 
elements in SNe la requires the WD to experience a phase of expansion to reduce densities 
in the outer layers and prevent the entire WD from burning to Ni. 



Intermediate mass elements such as Si, S, and Ca are produced in a wide range of 
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burning conditions (10^ g cm~^ > pbum > 4 x 10^ g cm~^ and 5 x 10^ K > T > 3 x lO^K). 
Similar abundance ratios for these elements are produced by both explosive O burning and 
incomplete Si burning. At the high end of this temperature and density range, silicon 
continues to burn after both carbon and oxygen are completely consumed, but the Si does 
not burn all the way to NSE. 

In regions with the lowest burning temperatures and densities [pbum < 4 x 10^ g cm~^ 
and T < 3 X lO^K), explosive carbon burning occurs generating ash that consists primarily 
O, Mg and Ne. The presence of these elements and the absence of C indicates that nuclear 
burning has occurred, since unburned regions of the SN will contain only C and O from the 
progenitor. Multiply-ionized atoms from heavier elements may be found in the outer layers 
from primordial sources. 



5.3. Location of the Photosphere and Line- forming Region 

We define the photosphere as the region within the expanding envelope where continuum 
photons are produced at a given time and wavelength. The photospheric radius is not 
constant, but is a function of wavelength. The specific location is a function of opacity in 
the envelope and the energy deposition rate from radioactive decay. As expansion reduces 
the optical depth in the outer layers, sequentially deeper layers of explosion products are 
revealed. A time series of spectra will trace the chemical composition from the outer regions 
toward the core as the photosphere recedes in mass space. NIR observations are particularly 
effective for recording the chemical structure of SNe la because the optical depth for most 
lines is less in the NIR than it is at s horter wavelengths. As a result, a greater radial distance 



can be probed with each spectrum (jWheeler et al.lll998l ). 



Absorption lines found in spectra from SNe la are formed in regions of expanding ma- 
terial that lie between the photosphere and the observer. These line-forming regions in SNe 
la are extended in radial space and hence in velocity space. To define the location of the 
line-forming region in radial space we use the Doppler shifts of the absorption components 
of individual spectral features to measure the expansion velocities. The radial extent of 
each line-forming region covers the full range of Doppler velocities between the absorption 
minimum and the detection limit in the blue wings of each absorption feature. 

During the first few weeks after the explosion, high temperatures within the envelope 
insure that the primary opacity source in the NIR is Thomson scattering with some contri- 
bution from free-free emission. In a rapidly expanding, scattering dominated envelope, all 
lines contain an emission component due to a non- vanishing line source function whether or 
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not they display an obvious P-Cygni profile. Emission photons from each line will reach the 
observer from different parts of the expanding envelope with different observed velocities. 
These emission photons contaminate the absorption feature and make Doppler measurements 
unreliable between the rest wavelength of the line and the wavelength of the absorption min- 
imum in each spectral feature. 

During this phase, while the photosphere and line-forming region are close, the mini- 
mum observed Doppler velocity is determined by the location of the photosphere. Because 
expansion velocity is proportional to radius, as the photosphere moves deeper with time, 
progressively slower material is exposed in the line-forming region. As a result, the mea- 
sured velocity at the absorption minimum for each line will diminish with time until there is 
insufficient abundance to produce a detectable signal. The measured range of velocities from 
the absorption minimum to the blue wing of the feature represents the physical location of 
the line-forming region in the ejecta. 

In the case where the line-forming region containing a particular ion does not extend 
below a certain radius, continued expansion will separate the photosphere from the line- 
producing layer. The observed Doppler velocity in all spectra obtained after separation 
occurs will be the velocity at the inner edge of this layer, and the measured Doppler velocity 
at the absorption minimum will no longer change. The absorption feature will subsequently 
decrease in strength as the column density drops due to the expansion of the line-forming 
region. Since the photosphere is detached from the absorbing region, lines arc formed far 
outside the last scattering radius for continuum photons and the photosphere appears to be 
nearly point-like. 

We note that during very early epochs (> 10 days before Vmax), the envelope is actually 
cooler than it will soon become, when it warms due to increased energy deposition from 
radioactive decay. The cooler material has a lower optical depth, causing the photosphere 
to temporarily recede far enough toward the center that some of the intermediate mass 
materials beneath the extreme outer layers are revealed. Doppler velocities are relatively 
low since the photosphere is closer to the center of mass and higher pressures in the interior 
near the photosphere will produce higher ionizations. Thus multiple ionizations stages from 
intermediate mass elements may be detected in a single spectrum due to the large range of 
radial distance probed by NIR spectra. 

Later than approximately ten days before V^ax, the expanding material continues to 
be heated by radioactive decay and the increasing temperature increases the opacity and 
moves the photosphere away from the center in both mass and velocity space. From this 
time forward, the photosphere begins the monotonic recession in mass and velocity space 
that characterizes spectral evolution during the next 15-20 days, as the data in this paper 
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demonstrate. 

A few days after Vmax , line-blanketing from thousands of Fe-group lines begins to dom- 
inate the spectrum at certain wavelength regions. The increase in total opacity moves the 
radius at which features are formed well above the receding photosphere formed by electron- 
scattering. A larger radius increases the effective area from which the observed flux is emit- 
ted. For wavelen gth regions where th is line-blanketing occurs, increased flux levels mimic 



emission features (jWheeler et al.lll998l ). There are thousands of these lines and they are not 
individually resolved but they effectively raise the continuum level through the wavelength 
regions that they occupy by creating a temporary effective photosphere at larger radii. Figure 
[2T] compares model data with four spectra from SN 2005am. The model resolution creates 
a series of peaks rather than a smooth profile but Doppler shifting will blend the individual 
lines and the models clearly predict the behavior observed in the data where a pseudo photo- 
sphere is created at extended radii that increases the continuum flux in wavelength regions: 
0.9-1.2, 1.5-1.85, and 2.05-2.45 fim. 

At 20-30 days after Vmax, the absorbing material is diluted by expansion to the point 
where it is no longer possible to define a photosphere. This is a transitional phase between the 
photospheric epoch and the true nebular phase. In some wavelength regions, line-blanketing 
from the thousands of Fe-group lines continues to distort the continuum. 

The Ca II IR triplet forms the strongest single feature in each of our spectra. Extremely 
low excitation values for these lines suggest that they are likely to create absorptions of 
this magnitude at very low Ca II abundances. Since solar values in circumstellar material 
would provide sufficient calcium to generate these features, it is not expected that the P- 
Cygni profile features found with absorption minima near 0.82 fim are produced by freshly 
synthesized Ca II. 



6. Sumniciry 

We present forty-one high quality NIR spectra from SNe la obtained between two weeks 
before and ten weeks after Vmax- We show that the Fourier transform is an excellent method 
for smoothing SN la data to allow more consistent measurements and better comparisons 
between spectra without compromising the important features. The procedures we describe 
for fitting continua and measuring absorption features provide consistency for measurement 
and analysis of a very large number of spectral features from spectra with a wide range of 
quality and resolution. We estimate relative line strengths for the strongest lines from ions 
expected to be found in spectra from SNe la in this wavelength range and provide tables of 
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the estimates as a resource for exploring NIR spectra. 

By comparing true time series of spectra with "snapshot" spectra from many individual 
SNe la, we show that the evolution of spectral features is consistent for NIR spectra from 
normal SNe la. Because geometrical dilution of the SN reveals deeper layers within the ex- 
panding envelope, our sample of forty-one spectra obtained between -14d and +75d provides 
an excellent record of the physical properties of SNe la from the outer layers to deep within 
the envelope. 

The spectra disclose that burning products are distributed in well-defined, radial layers. 
This is an important result that constrains the progenitor evolution and the burning physics. 
Table H] and Figure [19] show the measured Doppler velocities for O I, Mg II, Si II, S I, 
Ca II, and Fe II. Burning products from the lowest burning densities appear at the highest 
velocities (O I and Mg II at ~ 10,000-15,000 km s~^) with well-defined minima. High 
velocities correspond to the outer layers of SNe la and clearly-defined minimum velocities 
indicate that the line-forming region containing O and Mg has a distinct inner edge. Si II, 
S I, and Mn II are found at somewhat lower velocities (~ 8,000-12,000 km s~^) and they 
also have clear velocity boundaries to the line-forming region. The agreement in spectral 
features from a large number of SNe la eliminates the possibility that observations of finite 
line-forming regions are due to clumps of material. 

Ca II is only detected in the IR-triplet despite expectations that other Ca II lines would 
be present. Due to the very low excitation values for these lines, this result implies that 
detection is not due to freshly synthesized material. Fe II is detected at a wide range of 
velocities (~ 2,000-12,000 km s~^) with no apparent minimum. 

The distribution of burning products and radially stratified chemical structure is con- 
sistent with the explosion occurring in a C/0 WD with a radial density gradient and lower 
densities on the outside ( §5.2p . Our results agree with models that predict explosive carbon 
burning will occur in the extreme outer layers of SNe la. Oxygen burning and partial silicon 
burning appear to form intermediate layers and some Si must be completely burned to pro- 
duce Fe detected at extended radii. There are small overlaps between adjacent layers, but 
our data eliminate the possibility of large scale mixing of lower mass elements to the center 
or Fe-group elements to the extreme outer layers. 

C, O, and Mg are only detectable by spectral observations up to a few days after Vmax, 
so it is useful that there are twenty-two spectra in our sample obtained at -|-5d or earlier. 
We do not find carbon in the spectra at any epoch but we do find carbon burning products 
(O and Mg) occupying the same physical space in the outer layers of the SNe. There are 
strong hues from C I in the NIR that have produced detections in subluminous SNe la, but 
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C I is not found in our sample. (C II lines are expected to be weaker than C I and are also 
not detected.) O I and Mg II are both products of explosive carbon burning and they are 
detected in nearly the same velocity space in the outer layers. The line strengths for C and O 
will be similar for the same ionization stages, so if the O I in our spectra is du e to unburned 



material, then C I should also be present in detectable amounts and it is not (IMarion et al. 



20061). 



These results strongly suggest that the entire progenitor is consumed in normal SNe la 
which implies that nearly the same amount of fuel is used to power each explosion and helps 
to explain the uniformity in SNe la observations. These results constrain burning physics 
and models on a physical rather than phenomenological basis. 

The size of this sample also facilitates the search for secondary relationships with other 
data libraries that may further constrain SN la physics and improve the effectiveness of SNe 
la as cosmological distance indicators. Well-sampled light curves for the SNe in our sample 
will allow comparison of well-studied photometric parameters with NIR line strengths. NIR 
observables can now be compared to empirical patterns such as the relationship between peak 
brightness and light curve shape, B-V color differences, secondary peaks in the IR, and host 
galaxy morphology. This sample may enable accurate "K corrections" for photometric ob- 
servations in the J, H, and K bands. NIR spectra may also provide insights into areas of SNe 
la physics such as progenitor ZAMS mass, met allicity, rotation, asymmetry in the explosion, 
and gravitational influence on the detonation ( 



IQQsl : iHoflich et al.l[2002l : IYoou fc Langei1l2004al jbh. We hope that this sample will truly 



joflich et al.lll998l : IWang. Hoflich fc Wheeler 



open 



the door to use of NIR data for productive SNe la research. 
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A. Discussion of Individual Features and Line Identifications 

This section discusses line identifications for specific features in the spectra. Line iden- 
tifications and measured Doppler velocities are provided where possible. Individual lines are 
referred to by their rest wavelengths. The spectra are grouped by epoch using the same 
groups as in §11 

Velocities measurements discussed in this appendix can be found in Table [Hand Figure 

m 

A.l. A Very Early Spectrum (14 days before Vmax) 

Discussion in this section makes reference to Figure [TT] which shows the spectrum from 
SN 2002fk which was obtained fourteen days before Vmax- Features are labeled in the figure 
for which we have confident identifications. 

He I is not detected in this spectrum from either of the strongest He I lines in the NIR 
found at 1.0830 or 2.0581 /xm. Possible associations for He I include the absorption feature 
found near 1.05 /zm which corresponds to a Doppler velocity of 8,600 km s^^ for the 1.0830 
fim line. This is an unrealistically low velocity because if He I is present, it is expected to 
be found only in the extreme outer layers of the supernova at high velocities. The 1.05 fim 
feature has a plausible alternative identification as Mg II from the 1.0927 fim line and it is 
labeled as such in the figure. The 2.0581 fim He I line (which is expected to be 6 — 12 times 
weaker than the 1.083 fim line depending on temperature), can possibly be associated with 
either of two weak features at 2.00 and 2.04 fim with velocities of 9,700 or 2,500 km s~^ 
which are inconsistent with 8,600 km s~^ if the strongest He I line at 1.0830 produces 
the 1.05 /im feature. 

C I is predicted to generate strong lines at 1.0693, 0.9093, 0.9406, and 1.1754 /im in order of 
expected strength (see Tabled]). Additional discussion of C I is found in §4.31 and illustrated 
in Figure [201 Very weak features can be seen in this spectrum by optimistic observers at 
Doppler velocities of 10,000-12,000 km s~^, but no confident detections can be made and 
the velocities are not consistent within each spectrum. Some models for SNe la predict that 
C I will be present in the extreme outer layers, but we are unable to identify C I in this, or 
any other spectrum in our sample. 
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O I is likely from the 0.9264 /xm line contributing to the absorption feature found near 
0.89 fim at 11,500 km s'^ in a blend with Mg II (0.9227 /im). The 0.8446 O I line may 
also contribute to the absorption from the Ca II IR triplet near 0.81 but this wavelength 
is at the extreme edge of the spectrum and it is difficult to evaluate the presence of small 
features in this region due to the strength of the Ca II feature. The O I 0.7773 /im line is 
expected to be 3-20 times stronger than the 0.9264 line and the 0.7773 /xm line is detected 
in all early spectra from which we have coverage of the wavelength region down to 0.7 /im 
but this spectrum does not cover that range. 

Mg II produces strong features that are clearly detected in this spectrum (and labeled in 
Figure [TTl) from 0.9227 /im, found near 0.89 /im at 10,300 km s~^, and 1.0927 /im, found near 
1.05 /im at 11,300 km s~^. The measured velocity for the 0.9227 /im line is 1,000 km s~^ less 
than the measured velocity for the 1.0927 /im line due to intrusion on the blue side of the 
Mg II absorption from the emission component of the very strong Ca II IR triplet. This line 
is also likely to be blended with O I at 0.9264 /im. The high quality of this spectrum reveals 
a probable detection near 0.97 /im for the Mg II line at 1.0092 /im with a Doppler velocity 
of 10,900 km s~^. This line is predicted to be 5-20 times weaker than the 0.9227 /im line 
which is consistent with the observed feature strengths. The absorption feature near 1.60 
/im receives contributions from Mg II at 1.6787 /im with 11,600 km s^^ in addition to Si II 
(1.6930 /im) and Co II (many lines 1.57-1.77 /im). Weak evidence can be found for Mg II 
from the 0.9632 and 2.1369 /im lines at similar velocities. 

Si II from the strong line at 0.9413 /im is the likely source of the feature found near 0.91 
/im at a Doppler velocity of 11,300 km s~^. This velocity is similar to the velocities measured 
for Mg II which is consistent with models for this epoch. The 1.6930 Si II line may be part 
of the 1.6 /im absorption feature but the velocity to the minimum of the feature is nearly 
18,000 km s~^ for this line. This ion may also contribute to the feature found near 1.09 /xm 
from the 1.1311 /im line in a blend with Si III and Fe III. However, Si II at 1.1311 /im is a 
much weaker line than the 0.9413 /im line, and the velocity does not correlate with other 
Si II lines. If there is sufficient Si II to produce a feature of this magnitude, then detections 
should also be possible from lines at 1.1737 and 1.7183 /xm. There is no evidence in this 
spectrum for these other Si II lines. 

Si III can be present in the same spectrum with Si II if the photosphere is deep within the 
envelope so that temperature and pressure differences between absorbing regions near the 
photosphere and in the outer layers contain different ionization stages of the same element 
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(see §5.3p . That scenario suggests that Si III will be found in slower and hotter line- forming 
regions near to the photosphere. Due to a relatively low excitation minimum and high 
oscillator strength, the Si III 0.9324 fim line is predicted to be ~ 40 times stronger at 10,000 
K than the next strongest Si III line in the NIR. This line may contribute to the feature 
observed near 0.91 /zm at a Doppler velocity of 8,400 km s~^ in a blend with Si II and Mg II. 
A table of estimated line strengths using a temperature of 10,000 K shows change in the 
relative line strengths for Si III compared to 5,000 K (see Table [6]). For example: the fact 
that the Si III 1.0526 /im line is not detected in this spectrum is consistent with higher 
temperatures in the Si III line-forming region. Other possible identifications for Si III are 
the feature at 1.22 fim from the 1.2541 fim line (a blend of 1.2523 and 1.2601 fim) at 9,600 
km s~^ and 1.09 fim from 1.1341 fim at 10,500 km s~^. 

The Ca II IR triplet is strong in this spectrum from SN 2002fk with a Doppler velocity 
of 12,300 km s~^. The velocity measurement is based on a blend of the 0.8498 and 0.8542 fim 
Ca II lines, giving a mean wavelength of 0.8538 fim, and measured with respect to the bluer 
and deeper component of the absorption feature. Resolution is sufficient to separate this 
pair from the other member of the triplet found at 0.8662 /im. The shape of the absorption 
feature is consistent with estimates that the 0.8538 blend is about twice as strong as the 
0.8662 line. 

The very strong P-Cygni profile from the Ca II IR triplet may not be produced by 
freshly synthesized Ca ( §5.31) . There is no evidence in the spectrum for Ca II from the next 
strongest lines at 1.1876 or 0.9906 fim and this may indicate that Ca II produced in the 
explosion is not detectable. 

This SXD spectrum lacks coverage in the wavelength region to detect a possible high- 
velocity component from Ca II. 

Fe II lines are more difficult to compare based on our simple line strength estimates because 
the actual line strengths are affected by multiple transitions. As a consequence, the simple 
calculations in our tables for the relative estimated strength of each Fe II line may be less 
accurate than the estimates for smaller atoms. The Fe II 0.9998 fim line is a possible source 
for the feature found near 0.97 fim in this spectrum, but the velocity would be only 8,100 
km s~^. This velocity is similar to some of the conjectured Si III velocities. The 1.0500 ftm 
line from Fe II is not detected, but the 1.0863 line may be a part of the feature near 1.05 
ftm in a blend with Mg II. The 1.1126 ftm Fe II line is a possible contributor near 1.09 ftm 
in a blend with Si III. 
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Fe III at 0.9124 is possible in a blend with O I and Mg II but the Doppler velocity at the 
absorption minimum is very low at 7,100 km s^^. As noted in the discussion for Mg II in 
this section, this absorption feature is distorted by emission from the Ca II triplet and the 
true velocity may be higher. This ion has several lines from 1.60-1.67 fj^m that may add to 
the 1.6 /xm feature with the strongest line at 1.6697 fim having a velocity of 9,900 km s~^ 
at the minimum. The 1.1323 /im line may also form, or contribute to, the absorption near 
1.09 /zm at a velocity of 10,000 km s~^. 

Co II is detected at ~ 12,000 km s"^ from individual lines at 1.7772, 2.2205, 2.4596 , and 
2.3613 fim. The 1.7462 /zm line is not found in this spectrum. Numerous Co II lines between 
1.57-1.77 nm are likely contributors to the 1.60 /im feature. 

A. 2. The Photospheric Phase (Outer Layers) 
(-14d to +5d with respect to Vmax) 

During this epoch the photosphere is expected to be in the outer layers of the SN and 
the line-forming regions will be close to the photosphere. Discussion in this section makes 
reference to Figure [12] showing eighteen spectra obtained between -14 and +5 days with 
respect to to Vmax and also to Figure [13] that displays four representative spectra from this 
group: SNe 2004bw (LRS obtained at -9d), 2005am (SXD obtained at -4d), 2000dm (SXD 
obtained at +0d), and 2001bf (SXD obtained at +4d) with features labeled for which we 
have confident identifications. These figures only display the spectra in the wavelength region 
0.65 - 1.35 /xm to permit detailed inspection of the region where most of the features occur 
in this epoch. The complete spectra from 0.7-2.5 fim can be viewed in Figure [TOl 

C I is not detected from any of these spectra. The strongest C I lines are found at 1.0693, 
0.9087, 0.9406 and 0.9639 fim, listed in order of expected strength, (see Figure [20] and §4.31 
for additional discussion of C I in our sample). 

O I at 0.7773 forms the feature near 0.75 fim detected in all LRS spectra in this group (the 
SXD spectra do not cover the wavelength region of this line). Table [1] shows the Doppler 
velocity of this line is measured at 11,300-13,500 km s~^ except for velocities of 9,300 and 
9,400 measured in the two spectra from SN 2002cr. (Note that the spectra from SN 2002cr 
exhibit lower velocities than other spectra for all explosion products.) The O I feature is 
strong in all spectra from this group that include the appropriate wavelength region. 
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The 0.9263 /iin line from O I is expected to be 3-17 times weaker than the 0.7773 /im 
hne but it is hkely to contribute to the 0.89 /im feature in a blend with Mg II 0.9227 /im at 
velocities of 11,000-15,000. Possible evidence for the 1.1287 /im O I line (which should be 
half as strong as the 0.9263 /im line) appears near 1.09 /im in some of these spectra. 

The strongest features at 1.09 /^m occur in spectra from SNe 2002hw (-Id) and 2000dm 
(+ld) and correlate with the strongest absorptions at 0.89 /^m which is suggestive that O I 
contributes to both features. 

Mg II is found in this group producing strong features from lines at 0.9227 and 1.0927 /im 
and weaker but detectable features in higher S/N spectra from lines at 1.0092 and 1.6787 
/im. The Doppler velocities for the 1.0927 line are 11,100-14,500 km s~^ while the measured 
velocities for the 0.9227 line are ~ 1000 km s~^ slower as a result of the enormous P-Cygni 
emission from Ca II that impinges on the blue side of this absorption feature (Table S}. 
Stronger Mg II features are often correlated with higher measured Doppler velocities. As 
they were for O I, the lowest velocities are found in the spectrum of SN 2002cr. Figures [9] 
and [T2l show that there is an obvious correlation between strong spectral features from Mg II 
at 1.0927 /im and the strongest features from the 0.9227 /im line. The final spectrum (in 
temporal order) to show evidence of Mg II (SN 2000dk, +5d) clearly shows features from 
both of these strong lines. 

The Mg II line at 0.7890 /im may contribute to the 0.75 /im feature in a blend with O I 
but the velocities to the absorption minima of this feature would be 14,000-18,000 km s~^ 
and that is inconsistent with other Mg II measurements. The 0.7890 /im line may flatten 
the P-Cygni emission component of the O I line as in the spectrum from SN 2002cr. 

The 1.0092 /im Mg II line is expected to be 2 orders of magnitude weaker than 1.0927 
fim but it is detectable in most spectra in this group, forming a weak absorption near 0.97 
fim at velocities comparable to the other Mg II lines. This feature is detected in all but the 
noisiest spectra from this group and it often exhibits a double or triple bottom containing 
minima near 0.95, 0.97 and 0.99 /im with 0.97 /im being the strongest. The 0.95 and 0.99 /im 
features are smaller than the one sigma noise but are consistent from spectrum to spectrum 
in magnitude and location relative to the central feature. The central feature is probably 
Mg II at 1.0092 /im, while candidates for the others are Ca II at 0.9906 /im and Fe II at 
0.9998 /im. The fact that the feature near 0.97 /im becomes stronger after maximum light 
suggests that the primary source of absorption at this wavelength may be shifting from Mg II 
to Fe II. 

Velocity measurements of the Mg II line at 1.6787 /im and the 1.60 /zm absorption 
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feature correspond closely to the velocities of the other Mg II lines. This broad feature 
probably receives contributions from other ions, notably Si II and Co II. The absorption 
feature is particularly strong in spectra that exhibit other strong Mg II features, such as: 
SN 2005am (at -2d and +2d), SNe 2002hw (-Id) and 2000dm (+ld). That may indicate 
increased abundance of Mg II near max light. Mg II lines at 0.9632 and 2.1369 /im are also 
possible in weak features or blends. 

Inspection of the complete sample (Figure M> shows that all twenty-two spectra in our 
sample obtained at or before +2d show strong Mg II features near 0.89 and 1.05 /zm. Even 
spectroscopically unusual SNe la such as 2004bv exhibit distinct Mg II signatures. Measured 
velocities for Mg II are given in Table H] and plotted in Figure [191 

Si II from the 0.9413 /xm line is the probable source of the feature found near 0.91 fim. 
Note that this feature is distinct from the absorption found near 0.89 /im and it does not 
appear in all spectra. The measured velocities are between 10,200 and 14,100 km s~^. As 
noted in the discussion of the Mg II line at 0.9227 fim, impingement from the emission 
component of the Ca II triplet may affect the location of this feature, pushing it toward the 
red. The 0.91 /im feature, and thus the 0.9413 /xm Si II line, is not detected in spectra from 
SNe 2004bw (-9d), 2003W (-7d and -6d), 2002hw (-Id), 2000dm (+0d), or 2000dk (+5d). In 
general, the spectra without this Si II detection show very strong Mg II features that may 
obscure the Si II signal due to P-Cygni emission from Mg II at 0.9227 /im. 

The 1.6930 /im line from Si II is likely to contribute to the 1.60 /im feature but velocities 
to the absorption minima are 14,000-18,000 km s~^ which is higher than other measured Si II 
velocities including measurements during this epoch of the strong Si II features at optical 
wavelengths. The line at 0.7846 /im may contribute to the 0.75 /im feature in a blend with 
O I but velocities to the absorption minima are 12,200-16,400 km s~^ which is about 2000 
km s~^ higher than for the 0.9413 /im line. Possible evidence for the 1.1748 /zm line appears 
near 1.13 /zm in four of the spectra at plausible velocities of 10,100-11,900. However three 
of the four spectra with 1.13 /im features are the ones that do not have a 0.91 /im feature, 
making the Si II association unlikely. 

Si III is not detected in this group except for the earliest spectrum at -14d and discussed 
in the previous section. 

S I probably begins to contribute to the 0.89 /im feature after maximum light from 0.9223 
/im line at velocities of 10,000-12,000 km s~^ But the S I 1.0457 /im line, which is only 3 
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times weaker than the 0.9223 /xm hne, is not detected. 

S II is not detected in these data but all of the strongest hnes for this ion are on the blue 
side of 0.83 iim so their Doppler shifted positions will not be included in the SXD spectra 
that comprise the majority of spectra (and highest resolution) in this group. 

Ca II IR-triplet dominates the spectra in this group, forming the 0.82 /xm feature with 
velocities near 12,000 km for most SNe in our sample. The only exception is the spectrum 
from SN 2004bv (-4d) which may be a SN 1991T like object in which the Ca II is absent from 
early spectra but present, as in spectra from normal SNe la, at later times. The spectrum 
in this sample from SN 2004bv at +60d shows the Ca II and seems to be normal. 

The square bottom of the Ca II absorption region in these features is due to separately 
resolved contributions from two parts of the triplet. On the blue side of the absorption 

feature, we use a mean wavelength of 0.8538 /xm for the blend, weighted by oscillator strength, 
for the pair of lines at 0.8498 and 0.8542 fira. The 0.8662 /xm line can be easily identified in 
these spectra as a separate part of the feature at about one half the predicted strength for 
the 0.8538 fim line. To measure Doppler velocities, we use the location of the minimum at 
the blue side of the feature and compare it to the 0.8538 nm rest wavelength of the blend. 

The shape of the emission region forming the P-Cygni profile apparently refiects the 
shape of the absorption region. Flat or round shapes in the absorption component correspond 
to similar but inverted profiles in the emission component. 

High Velocity Ca II The spectrum of SN 2004bw (-9d) shows two distinct minima for 
Ca II with one at 13,400 km s^^ and a high-velocity component at 24,000 km s^^. Both 
spectra from SN 2003W (-7d and -6d) have Ca II velocities above 22,100 km s~^ without a 
lower velocity component. The spectra from SN 2002cr (-7d and -5d) appear to have a fading 
high velocity Ca II component at 18,800 km and a lower velocity Ca II absorption at 
9,700 km s~^. (Both spectra from SN 2002cr exhibit lower velocities for all measured features 
than other spectra at similar epochs.) SN 2000dn (-6d) also has a weaker Ca II feature near 
22,000 km s~^ and a strong Ca II absorption near 10,400 km s~^. All LRS spectra obtained 
at -5d or earlier appear to have high velocity Ca II components with velocities measured 
between 18,800-24,000 km s~^. Since the high- velocity features have rounded bottoms, we 
determine the Doppler velocities for the high-velocity Ca II feature using 0.8579 /xm as the 
mean wavelength for a blend of all three strong lines in the Ca II IR-triplet. 
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Other Ca II may be detected in some spectra near 1.14 /im from the 1.1876 /xm hne at 
velocities of 12,000-16,300 km s~^. The spectra in which this feature is the most obvious also 
have the flattest tops to the IR-triplet emission and that could be due to the next strongest 
Ca II line at 0.8921 fim at similar velocities. The feature attributed to Mg II near 0.97 nm 
may have a contribution from Ca II at 0.9906 /zm but the velocity would be ~ 8,000 km s~^ 
which is low. The weak absorption feature near 0.97 /im often exhibits a double or triple 
bottom with minima near 0.95, 0.97 and 0.99 /im with 0.97 being the strongest. The 0.95 
and 0.99 fim features are smaller than the one sigma noise but are present in most spectra. 
An absorption near 0.95 fim created by Ca II at 0.9906 fim would have a Doppler velocity 
of 12,300 km s'^ 

Mn II produces no unambiguous detections but the Mn II lines at 0.7414 and 0.7342 /im 
are a possible source for the feature found near 0.70 /im in LRS spectra. This identification 
would require velocities of 14,000-17,000 km s~^ which are higher than those found for most 
of the other species. The strong Mn II line at 0.9406 /im may contribute to the feature near 
0.91 /im either alone or in a blend with Si II at velocities that are ~ 300 km s~^ less than 
values given for Si II at 0.9413 in Table HI The 0.8410 line is also possible in the blend 
dominated by the Ca II IR-triplet. 

The steep increase in flux found near 1.55 /im is used by some models to define the 
transition from partial to complete silicon burning. This feature is suggested in the spectra 
from this group as early as +0d but it becomes distinctive by -|-4d. 

Fe II at 0.9998 /im probably contributes to the 0.97 /im absorption that we have associated 
with Mg II (1.0092 /im). The fact that this feature becomes stronger after max light suggests 
that the primary source for this absorption is shifting from Mg II to Fe II. Doppler velocities 
for the Fe II association with this feature begin above 10,000 at -4d and follows as the 
photosphere recedes into regions with lower expansion velocities to about 7,000 km s^^ +2d. 

The velocities required to match the 0.7712 /im Fe II line with the 0.75 /im absorption 
feature in these spectra are inconsistent with Fe II velocities measured for the 0.9998 /xm 
line and with other species found in this group. There is no evidence of Fe II from lines at 
0.7462 or 1.0500 /im. The 1.0863 /im line may contribute to the 1.05 /im feature in a blend 
with Mg II but this line is expected to be ~ 20 times weaker than some of the other Fe II 
lines that are not detected. 

Fe-group lines are complicated by many potential transitions at the same wavelength 
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and nearly the same strength within each multiplet so the simple formula we have been 
using to estimate line strength may not be accurate. In this group, the 0.9998 /im line is the 
only Fe II line detected, although the calculated line strengths suggest that lines at 0.7462, 
0.7712, and 1.0500 /im should also be detected. 

Co II probably contributes to the 1.60 /zm feature with numerous Co II lines between 
1.57-1.77 fim. In a few spectra (SNe 2002hw (-Id), 2000dm (+ld), and 2005am (+2d)) the 
strongest Co II line (1.7772 /im) may be separately resolved at about 11,000 km s~^. 

Fe-group blends begin to affect the spectra in this group as line-blanketing from thou- 
sands of Fe-Group lines results in pseudo-emission features found near 1.52 /im and between 
1.6-2.1 /im. Model predictions for this effect are compared to the multiple spectra obtained 
from SN 2005am in Figure [21] and discussed in Section 15.31 Broad humps, that we attribute 
to unresolved Fe-group lines raising the continuum, are also observed in the spectra from 
1.6-2.1 /im and from 2.1-2.3 /im. The bump near 1.25 /im is very strong early in this group 
(possibly Si III), weaker through the middle and strong again in the later spectra. These 
features appear in all spectra from this group. 

A. 3. The Extended Photospheric Phase (Intermediate Layers) 
(-|-4d to -|-24d with respect to Vmax) 

During the phase represented by this group, the shape of the continuum becomes influ- 
enced by line-blanketing from Fe-group lines which extends the effective photospheric radius 
(see §5.3p . Discussion in this section refers to Figure [T3] showing ten spectra from SNe la 
obtained between +4d and +24d with respect to to Vmax and also to Figure [T5l that displays 
four representative spectra from SNe 2000dk (+5d), 2001bg (+10d), 2004ab (+18d), and 
2004da (+24d) with the features labeled for which we have confident identifications. 

O I produces the strong absorption near 0.75 /xm in the LRS spectrum of SN 2002ha (+8d) 
and this line may also be responsible for the small features at this wavelength up to +18d. 
The 0.9263 /im line from O I is expected to be 3-17 times weaker than the 0.7773 /im line but 
it is likely to be a part of the 0.89 /im feature in a blend with S I at least early in the epoch 
represented by this group. The measured Doppler velocities for the strongest detections of 
O I are all 9,800 km s~^ or higher. 

Alternative explanations for the 0.75 /xm feature are Fe II or Mn II, but the Doppler 
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velocities for these lines do not fit velocities measured for other possible lines from these 
ions. 

Mg II is found in this only in the spectrum of SN 2000dk at +5d which is the oldest 
spectrum int the sample to show strong Mg II absorption. Mg II is the primary contributor 
forming the distinctive notch near 1.60 /xm. This spectrum also has strong Mg II absorption 
from the 0.9227 and 1.0927 fim lines. The 1.0092 /xm Mg II line is expected to be 2 orders of 
magnitude weaker than 1.0927 /xm but at +5d or earlier it may contribute to the absorption 
near 0.97 /xm either alone or in a blend with Fe II. 

Si I is not detected in our data although this ion produces many strong lines in the NIR 
including lines at 1.2032, 1.0827, and 1.0585 fim. These lines may help form the features 
found near 1.03 and 1.06 /xm in the oldest spectrum from this group (SN 2004da at +24d). 

Si II is possible with the 0.9413 fim line contributing to the 0.92 /xm absorption found in 
some spectra in this group. The Si II 1.6930 fim line may also be part of the 1.60 fim feature. 
The 1.1737 fim line from Si II is not detected so it is unlikely that the 1.1311 fim line is 
responsible for the feature found near 1.12 fim in some of the spectra. The 1.12 fim feature 
does not occur in the same spectra as the 0.92 fim feature which eliminates the possible 
association of the 1.12 //m feature with Si II. 

SI at 0.9223 fim contributes to the 0.90 fim feature in all spectra in this group and probably 
replaces Mg II as the primary source of the absorption after +5d post-maximum. Velocities 
for this ion begin at 10,300-13,700 km s"^ at +0d and diminish to 8,900 at +18d. After 
+10 days, some spectra do not have any absorption feature at this location. The 1.3809 
fjim S I hue is too far into the region of high atmospheric opacity to be detected. Other 
possible associations can be made for S I at 0.9423 and 1.1349 fim, but they both have 
serious questions. The feature found near 1.10 fim is possibly associated with the 1.1349 
fim line and the velocity shift is plausible of 8,500 at +5d and 7,000 at +18d. However this 
feature is near an order transition in SXD spectra which casts suspicion on any features at 
this location from SXD spectra (in particular the spectrum from SNe 2000dk (+5d)). Since 
LRS spectra do not have order transitions, the weak feature near 1.1 fj,m is more likely to be 
real in the LRS spectra. The S I 0.9423 fxm line may also produce features in a few spectra, 
but the 0.92 fim absorption feature does not appear in the same spectra that produce other 
probable S I identifications. 
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Ca II IR triplet velocities remain well above 10,000 km for the spectra in this group. 
Separate features can be detected in the spectra from this group for both the 0.8358 /xm blend 
and the 0.8662 /xm line as described in the previous section. The shape of the absorption 
feature is consistent with estimates that the 0.8538 fim blend is about twice as strong as the 
0.8662 /im line. 

Other Ca II is not detected from the lines at 1.1876 or 0.9906 /im or from other Ca II 
line. The weak features found in some spectra near 1.10 fim are possibly associated with 
the 1.1876 /im line but the velocity would be more than 21,000 km s~^ which is inconsistent 
with other photospheric velocities during this phase. That makes it unlikely that the Ca II 
line at 0.9931 /im, which is estimated to be ~ 20 times weaker, provides a contribution to 
the absorption found near 0.97 /im. 

Mn II at 0.9447 /im is a possible detection in the spectra from this group that contain 
features near 0.92/im. This feature may be attributed to either Mn II or Si II alone or in a 
blend of both. A rest wavelength of 0.9447 /im would produce velocities of 8,500-12,500 km 
s~^ for Mn II, or about 1000 km s~^ greater than velocities listed for Si II at 0.9413 /im in 
Table HI However Mn II produces more than 90 strong lines between 0.9309-0.9509 /xm so 
it is not possible to measure a definitive velocity. These lines do not appear in Tables [5] or [6] 
because they are slightly weaker than the 0.9447, 0.9388, and 0.9486 /xm lines that do appear 
in the tables. The strongest Mn II lines at 0.7414 and 0.7342 /im may form the features 
found with a similar velocity separation near 0.72 /im, but there are only two spectra with 
which to compare these lines. An alternative identification for those features is Fe II. 

The steep increase in flux found near 1.55 /im may define the transition from partial 
to complete silicon burning. This feature is present in all spectra from this group but it is 
less well defined in the spectra from SNe 2005am (+14d) and 2004da (+24d). This does not 
appear to be an evolutionary effect because the features is present in many later spectra. 

Fe II at 0.9998 /im is the likely source of the absorption found near 0.97 /xm in all spectra 
from this group through +18d. Before +5d (see the previous section) this feature may receive 
a contribution from Mg II at 1.0092 /im, but the absorption is strongest from +5d to +10d 
when it is likely to be entirely due to Fe II. Velocities for Fe II measured from the 0.9998 
/xm line are near 9,000 km s~^ at +0d and diminish to less than 3,000 by +18d. Fe II is also 
likely to be the source of absorptions near 0.75 and 0.72 /im after +10d from the two strong 
lines at 0.7712 and 0.7462 /xm at about 8,500 km s"^ in the spectra from SNe 2004bk (+12) 
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and 2002ef (+16). The Fe II line from 0.9998 /xm appears to diminish in strength after about 
+10d but it becomes strong again about one month post-maximum. 



Co II forms distinct features in the later spectra in this group. Spectra obtained at +10d 
or later show distinct absorption at the location of many Co II lines, most prominently: 
1.7772, 1.7462, 2.2205, 2.4596, and 2.3613 /im at velocities near 9,000 km s"^ (see Figure 
[T7I) . Possible identification of Co II lines from 1.7772 and 1.7462 fim can be made in the 
spectrum from SN 2000dk as early as +5d. 



A. 4. Transitional Phase (Inner Layers) 
(+37d to +75d with respect to Vmax) 

This epoch is a transitional phase between the photospheric era and a true nebular 
phase. Discussion in this section makes reference to Figure [16] showing seven spectra from 
SNe la obtained between +37d and +75d with respect to Vmax and Figure [17] which includes 
13 spectra from our sample and two additional spectra from other sources. In Figure [T7] 
the positions of strong Co II and Fe II lines are marked so that they correspond to 11,000 
km at the top spectrum and 6,000 km s~^ at the bottom spectrum. This reflects to the 
reduction in measured velocities for these spectra . The data plotted in yellow (SNe 1998aq 



and 1998bu) were obtained by iRudy &: Puetterl and have not been previously published 



Co II lines are easily detected in absorption as early as two weeks post-maximum and Fe II 
lines (in addition to the 0.9998 fim which appears as early as +0d) begin to appear about 
one month post-maximum (see §3.1.31 and 13.1.41) . 

This group displays a lot of new structure relative to earlier NIR spectra and the features 
are remarkably similar from one spectrum to the next. The spectrum of SN 2004bv at +60d 
fits well with the others in this group although the spectrum from this SN at -5d was 
unusual (see Figure [T2l) . This is the same time frame in which unreddened SNe la have a 
very low dispersion for B-V colors (Lira, 1995). The common assumption about this era is 
that expansion has reduced the opacity beyond the point where a distinct photosphere can 
be defined and indiv idual absorption features can no longer be produced. Recent work by 



Branch, et al.l (l2008l ) suggests spectra from this epoch can be modeled by absorption features 



from permitted lines. 



Ca II IR triplet velocities remain near 10,000 km s ^ which is consistent with this feature 
being produced by Ca II that was not synthesized in the explosion (see §5.3p . Separate 
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features can be detected in many spectra from this group for both the 0.8358 /im blend and 
the 0.8662 fim line. 

The steep increase in flux found near 1.55 /im may define the transition from partial 
to complete silicon burning. This feature is present is all spectra from this group but is more 
sharply defined in some spectra compared to others. As shown in Figure [161 the spectrum 
from SN 2004bv at +60 displays this feature but it is smeared out in the spectrum of SN 
2001fe, also obtained at +60d. 

Fe II lines at 0.7712, 0.9998, and 1.0863 /xm appear in absorption beginning about 30 days 
after maximum light and in all spectra in this group (Figure [T71) . The Fe II 0.9998 /im line 
is present is spectra from about +0d to +18d, it is not detected from +18d to +24d, and 
then reappears after +30d (see Figure [TOj) . No other detections from individual Fe II lines 
are reliable from NIR spectra prior to this later phase. 

Co II lines are clearly detected in wavelength regions where a pseudo photosphere is created 
at extended radii (1.5-1.85 and 2.05-2.45 /im; see §5.3p and Figure [T71 Strong absorption 
features are detected from many Co II lines. The LRS spectra from SNe 2003cd (+45) and 
2004bv (+60) as shown in Figures [16] only marginally resolve Co II absorptions for 1.7772, 
1.7462, 2.2205, 2.4596, and 2.3613 lines and the group of Co II lines from 1.57-1.67 are not 
resolved at all. Figure [17] shows how the locations of these lines shift to longer wavelengths 
with increasing time of observation relative to maximum light. The positions of strong Co II 
and Fe II lines are drawn on the figure so that the indicated velocities are 11,000 km s~^ at 
the top spectrum and 6,000 km s~^ at the bottom spectrum. 

A. 5. Summary of Line Identification by Ion for the Entire Sample 

Velocity references are compiled in Tabled] and Figure UM 

He I is not detected from either of the strong lines at 1.0830 or 2.0581 /im. Unrealistically 
low Doppler velocities are required to align the 1.0830 fim line with the absorption feature 
found near 1.05 fim. The 2.0581 line (which is expected to be 6-12 times weaker than the 
1.083 line depending on temperature), is not observed. 
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C I is not detected. Speculative identifications can be made for lines at 0.9406, 1.1331, 
1.1756 and 1.6890 /xm in the very early spectra, but if these lines produce absorption features, 
the much stronger lines at 1.0693 and 0.9087 should also be detected and they are not. At 
each location of possible C I detection there are plausible alternative identifications (see 
Figure [20]). 

O I forms a strong feature from the 0.7773 fim line producing an absorption near 0.75 /xm. 
This feature is found in all LRS spectra up to +10d and possibly as late as +18d (SXD spectra 
lack the wavelength coverage to include this feature). O I is also likely as a contributor to 
the 0.89 feature from the line at 0.9264 /zm. Doppler velocities are 9,300-15,700 km s~^ with 
most measurements near 11,000 km s~^. 

C II &; O II are not detected. All C II lines in the NIR are factors of 10^ — 10^ weaker 
than the C II at 0.6580 /im line which has been tentatively identified in a few optical spectra. 

Mg II is found in all spectra between -14d and +5d, producing strong features from lines 
at 0.9227 and 1.0927 /im. Velocities are 9,800-14,700 km s~^ with a few anomalies near 
9,000 km s~^. Measured velocities for the 0.9227 line are ~ 1000 km s~^ slower than for the 
1.0927 /im line as a result of the enormous P-Cygni emission from Ca II impinging on the 
absorption feature. Many spectra also contain features from weaker Mg II lines at 1.0092 
and 1.6787 /xm. 

Si I is not detected although some spectra exhibit weak features that could possibly be 
attributed to this ion from lines at 1.2032, 1.0827, or 1.0585 fj,m. 

Si II from the 0.9413 /im line produces the feature near 0.91 /im at velocities from 10,200- 
14,100 km s~^ between -14d and +5d. After +5d the feature appears intermittently as late 
as +18d, but this may be a blend with Mn II or entirely due to Mn II. The Si II 1.6930 
line contributes to the 1.60 feature. The 1.1737 line is not detected so it is unlikely that the 
1.1311 line is responsible for the feature found near 1.12 in some of the spectra. The 1.12 
feature does not occur in the same spectra as the 0.91 feature so the Si II identification is 
unlikely. 

Si III is possible in the earliest spectra from the 0.9324 fim line at the very low velocity of 
8,500 km s~^. Other possible identifications for this ion are the feature at 1.22 (from 1.2541, 
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a blend of 1.2523 and 1.2601, at 9,600 km s"^) and 1.09 (from 1.1341 at 10,500 km s^^). 
Si III can be present in the same spectrum with Si II if the photosphere is deep within the 
envelope so that temperature and pressure differences between absorbing regions near the 
photosphere and in the outer layers contain different ionizations of the same element ( §5.3p . 

Si at 0.9223 fim replaces Mg II as the primary source of the 0.90 feature after maximum 
light. Velocities begin at 10,300-13,700 km s"^ at +0d and diminish to 8,900 at +18d. After 
+ 10 days, there are a few spectra with no feature at this location. Other possible associations 
can be made for S I at 1.3809, 0.9423, and 1.1349, but they all have serious questions. 

Ca II IR triplet is the strongest feature in most spectra. Velocities are greater than 
10,000 km s~^ for most spectra, although there are exceptions such as SN 2001bf (+4d) for 
which the velocity is 7,600 km s~^ and this feature is very weak in the spectrum from SN 
2004bv at -5d. SN 200 Ibf has very low velocities for all identified lines and SN 2004bv may 
be 91T-like. Distinct individual features can be detected in most of the SXD spectra for both 
the 0.8358 blend and the 0.8662 line. The shape of the absorption feature is consistent with 
estimates that the 0.8538 blend is ~ twice as strong as the 0.8662 line. Very low excitation 
values for these lines suggest that absorption features from the Ca II IR triplet may not 
come from freshly synthesized material. 

High Velocity Ca II All LRS spectra obtained at -5d or earlier appear to have high 
velocity Ca II components with velocities measured between 18,800-24,000 km s~^. 

Other Ca II is not detected from the next strongest lines at 1.1876 or 0.9906 /xm or from 
other Ca II lines in any spectrum in our sample. Ca II appears to be absent from the freshly 
synthesized explosion products. 

Mn II at 0.9447 fim is possibly detected after maximum light. The strongest Mn II lines 
at 0.7414 and 0.7342 fim form features at similar velocities near 0.72 /im, but there are only 
two spectra with which to compare this line. 

Fe II is detected in post-maximum spectra from the 0.9998 /xm line at velocities near 9,500 
km s~^ at +0d, diminishing to less than 3,000 by +18d. This line may contribute to the 
feature in pre-maximum spectra attributed to the Mg II line at 1.0092 /im. Fe II is likely 
to be the source of absorptions near 0.75 and 0.72 after +10d. Fe lines are more difficult to 
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compare based on our simple line strength estimates because the actual line strengths are 
affected by multiple transitions. As a consequence, the simple calculations in our tables for 
the relative estimated strength of each Fe II line may be less accurate than the estimates for 
smaller atoms. In spectra from our sample obtained before +37d, the 0.9998 fim line is the 
only Fe II line detected, although the calculated line strengths suggest that lines at 0.7462, 
0.7712, 1.0500, and 1.0863 /xm lines should also be detected. 

Co II may be detected soon after maximum light in some of the high resolution spectra. 
After two weeks post maximum, Co II lines are clearly detected at velocities beginning near 
10,000 km s~^ and diminishing to 6,000 km s"^. Co II is found in wavelength regions where a 
pseudo photosphere is created at extended radii (1.5-1.85 and 2.05-2.45 jim). Co II is most 
obvious in the higher resolution (SXD) spectra and individual Co II lines are marginally 
resolved in many other later time spectra. 
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Table 1. Comparing the Effects of Fourier Transform Smoothing Parameters 



Smoothing 
Parameter 


Minimum (^m) 
Xrest = 1.0927 fim 


Velocity 
km 


Line Depth 
normalized 


Minimum (fim) 
Xrest = 0.9227 fim 


Velocity 
km 


Line Depth 
normalized 


39 


1.0463 


12,700 


1.00 


0.8859 


11,900 


1.00 


32 


1.0467 


12,500 


0.99 


0.8864 


11,800 


0.99 


27 


1.0474 


12,400 


0.97 


0.8866 


11,700 


0.96 


22 


1.0479 


12,300 


0.95 


0.8873 


11,500 


0.82 


16 


1.0479 


12,300 


0.87 


0.8894 


10,800 


0.86 


11 


1.0471 


12,500 


0.78 


0.8921 


9,900 


0.82 



Note. — 



These data are plotted in Figure |4] and discussed in i|2.1l 



Table 2. List of Spectra by Estimated Epoch Relative to 



SN 


Est. Epoch 


Obs. 




S/N 


S/N 


S/N 


S/N 


Rodshift 


Name 


wrt V^^^ 

vv ± u V max 


Date (UT) 


Res. 


0.95 — 1.10/xm 


1.10 — 1.30/im 


1.45 - 1.70/xm 


1.95 — 2.35/im 


Host (z) 


2002fk 


-14d 


Sep 19.4 


1200 


74 


67 


47 


25 


0.007125 


2004bw 


-9d 


May 30.4 


200 


74 


54 


17 


12 


0.021408 


2003W 


-8d 


Feb 2.4 


200 


104 


29 


47 


22 


0.020071 


2003W 


-7d 


Feb 3.4 


200 


144 


132 


41 


25 


0.020071 


2002cr 


-Id 


May 8.4 


200 


153 


133 


77 


40 


0.009590 


2000dn 


-6d 


Oct 1.5 


120 


37 


34 


7 


3 


0.032019 


2004bv 


-5d 


May 30.5 


1200 


91 


69 


54 


37 


0.010544 


2002cr 


-5d 


May 10.3 


200 


53 


37 


11 


4 


0.009590 


2005am 


-4d 


Mar 5.3 


1200 


20 


21 


18 


12 


0.007889 


2001en 


-4d 


Oct 8.4 


1200 


15 


19 


10 


7 


0.015871 


2002el 


-4d 


Aug 20.5 


750 


6 


7 


4 


4 


0.023300 


2001br 


-2d 


May 22.6 


750 


13 


18 


9 


7 


0.020628 


2004bl 


-2d 


May 8.3 


1200 


7 


7 


4 


3 


0.017319 


2001dl 


-Id 


Aug 12.5 


750 


9 


12 


7 


5 


0.020694 


2002hw 


-Id 


Nov 14.3 


1200 


12 


12 


6 


4 


0.017535 


2005am 


+0d 


Mar 9.4 


1200 


30 


32 


22 


16 


0.007889 


2000dm 


+0d 


Oct 1.2 


750 


15 


21 


10 


9 


0.015304 


2003du 


+ld 


May 9.3 


1200 


47 


31 


24 


15 


0.016764 


2001bf 


-|-4d 


May 21.5 


1200 


50 


51 


36 


29 


0.015501 


2000do 


+4d 


Oct 2.2 


750 


15 


10 


7 


6 


0.010864 


2004da 


+Ad 


Jul 8.5 


750 


8 


8 


6 


5 


0.016308 


2000dk 


+5d 


Oct 1.5 


750 


19 


14 


11 


7 


0.017439 


2005am 


+7d 


Mar 16.4 


750 


3 


2 


2 


3 


0.007889 


2002ha 


+8d 


Nov 14.2 


200 


60 


23 


13 


10 


0.014046 


2001bg 


-MOd 


May 22.2 


1200 


18 


10 


9 


7 


0.007122 


2002ef 


-l-lld 


Aug 20.6 


200 


22 


7 


5 


2 


0.023977 


2005am 


-|-14d 


Mar 23.3 


1200 


26 


15 


19 


13 


0.007889 



Table 2 — Continued 



SN 


Est. Epoch 


Obs. 


Avg. 


S/N 


S/N 


S/N 


S/N 


Rcdshift 


Name 


wrt Vmax 


Date (UT) 


Res. 


0.95- l.lOAtm 


1.10- 1.30/iin 


1.45- 1.70/iin 


1.95 - 2.35Atm 


Host (z) 


2004da 


+ 14d 


Jul 18.4 


750 


4 


2 


5 


5 


0.016308 


2004ab 


+18d 


Mar 7.5 


750 


41 


27 


45 


36 


0.005807 


2001en 


+18d 


Oct 30.3 


200 


45 


21 


18 


10 


0.015871 


2004da 


+18d 


Jul 22.5 


750 


35 


15 


30 


23 


0.016308 


2004bk 


+19d 


May 8.5 


200 


12 


8 


2 





0.023036 


2UU4aa 


+24d 


T 1 no A 

Jul 28.4 


750 


1 O 

18 


10 


15 


10 


0.016308 


2004E 


+35d 


Feb 21.5 


200 


42 


13 


6 


1 


0.029807 


2003cg 


+37d 


May 9.2 


200 


66 


56 


80 


47 


0.004130 


2001gc 


+40d 


Jan 14.5 


200 


28 


12 


8 


6 


0.019290 


2002fk 


+42d 


Nov 14.4 


750 


37 


13 


15 


10 


0.007125 


2004ca 


+46d 


Jul 28.6 


200 


56 


24 


13 


5 


0.017806 


2004bv 


+53d 


Jul 28.3 


200 


63 


45 


38 


16 


0.010544 


2001fc 


+60d 


Jan 14.6 


120 


45 


15 


12 


6 


0.013539 


2003du 


+75d 


Jul 22.3 


200 


41 


7 


8 


2 


0.016764 



Table 3. List of Spectra in Order of Discovery 



SN 


Est. Epoch 


Obs. 


Avt? 


S/N 


S/N 


S/N 


S/N 


Redshift 


Name 


wrt V^«^ 

vv ± u V max 


Date (UT) 


Res. 


0.95 — l.lO/xm 


1.10 — 1.30/ini 


1.45 - 1.70/xm 


1.95 — 2.35/im 


Host (z) 


2000dk 


+5d 


Oct 1.5 


750 


19 


14 


11 


7 


0.017439 


2000dm 


+0d 


Oct 1.2 


750 


15 


21 


10 


9 


0.015304 


2000dn 


-6d 


Oct 1.5 


120 


37 


34 


7 


3 


0.032019 


2000do 


+4d 


Oct 2.2 


750 


15 


10 


7 


6 


0.010864 


2001bf 


+4d 


May 21.5 


1200 


50 


51 


36 


29 


0.015501 


2001bg 


+ 10d 


May 22.2 


1200 


18 


10 


9 


7 


0.007122 


2001br 


-2d 


May 22.6 


750 


13 


18 


9 


7 


0.020628 


2001dl 


-Id 


Aug 12.5 


750 


9 


12 


7 


5 


0.020694 


2001cn 


-4d 


Oct 8.4 


1200 


15 


19 


10 


7 


0.015871 


2001en 


+18d 


Oct 30.3 


200 


45 


21 


18 


10 


0.015871 


2001fe 


+60d 


Jan 14.6 


120 


45 


15 


12 


6 


0.013539 


2001gc 


+40d 


Jan 14.5 


200 


28 


12 


8 


6 


0.019290 


2002cr 


-7d 


May 8.4 


200 


153 


133 


77 


40 


0.009590 


2002cr 


-5d 


May 10.3 


200 


53 


37 


11 


4 


0.009590 


2002ef 


+lld 


Aug 20.6 


200 


22 


7 


5 


2 


0.023977 


2002el 


-4d 


Aug 20.5 


750 


6 


7 


4 


4 


0.023300 


2002fk 


-14d 


Sep 19.4 


1200 


74 


67 


47 


25 


0.007125 


2002fk 


+42d 


Nov 14.4 


750 


37 


13 


15 


10 


0.007125 


2002ha 


+8d 


Nov 14.2 


200 


60 


23 


13 


10 


0.014046 


2002hw 


-Id 


Nov 14.3 


1200 


12 


12 


6 


4 


0.017535 


2003W 


-8d 


Feb 2.4 


200 


104 


29 


47 


22 


0.020071 


2003W 


-7d 


Feb 3.4 


200 


144 


132 


41 


25 


0.020071 


2003cg 


+37d 


May 9.2 


200 


66 


56 


80 


47 


0.004130 


2003du 


+ld 


May 9.3 


1200 


47 


31 


24 


15 


0.016764 


2003du 


+75d 


Jul 22.3 


200 


41 


7 


8 


2 


0.016764 


2004E 


+35d 


Feb 21.5 


200 


42 


13 


6 


1 


0.029807 


2004ab 


+18d 


Mar 7.5 


750 


41 


27 


45 


36 


0.005807 



Table 3 — Continued 



SN 


Est. Epoch 


Obs. 


Avg. 


S/N 


S/N 


S/N 


S/N 


Rcdsliift 


Name 


wrt Vmax 


Date (UT) 


Res. 


0.95- 1.10/xm 


1.10- 1.30Atm 


1.45 - 1.70Atm 


1.95- 2.35/xm 


Host (z) 


2004bk 


+19d 


May 8.5 


200 


15 


8 


2 





0.023036 


2004bl 


-2d 


May 8.3 


1200 


7 


7 


4 


3 


0.017319 


2004bw 


-9d 


May 30.4 


200 


74 


54 


17 


12 


0.021408 


2004bv 


-5d 


May 30.5 


1200 


91 


69 


54 


37 


0.010544 


2004bv 


+53d 


Jul 28.3 


200 


63 


45 


38 


16 


0.010544 


2UU4ca 


-|-4od 


Jul 28.0 


200 


56 


24 


13 


5 


0.017806 


2004da 


+4d 


Jul 8.5 


750 


8 


8 


6 


5 


0.016308 


2004da 


+14d 


Jul 18.4 


750 


4 


2 


5 


5 


0.016308 


2004da 


+18d 


Jul 22.5 


750 


35 


15 


30 


23 


0.016308 


2004da 


+24d 


Jul 28.4 


750 


18 


10 


15 


10 


0.016308 


2005am 


-4d 


Mar 5.3 


1200 


20 


21 


18 


12 


0.007889 


2005am 


+0d 


Mar 9.4 


1200 


30 


32 


22 


16 


0.007889 


2005am 


+7d 


Mar 16.4 


750 


3 


2 


2 


3 


0.007889 


2005am 


+14d 


Mar 23.3 


1200 


26 


15 


19 


13 


0.007889 
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Table 4. Measured Doppler Velocities for Identified Lines - Blueshift (km s ^) 



Name 


Epoch 


I 


I 


Mg II 


Mg II 


Mg II 


Si II 


S I 


Ca II 


Fc II 




^max 


0.7773 


0.9264 


0.9227 


1.0092 


1.0927 


0.9413 


0.9223 


0.8538 


0.9998 








1 1 1 nn 


QQOO 


1 osoo 


1 1 "^00 


1 osoo 




1 9fi00 




iUU^ u w 


Q 




1 '^000 


1 1 SOO 


1 "^400 








1 "^OOO 


1 OfiOO 


^UUO vv 


s 
o 




1 1 non 


QSOO 


1 9fion 










Q7nn 


^uuo vv 


7 


1 oouu 


1 1 nnn 

1 lUUU 


QSOO 


i4sno 

±4:OUU 


1 1 nnn 








1 9000 




7 


e/OUU 




QSOO 




1 n nn 


1 1 fino 




Qsnn 


Q7nn 


9nnnrin 


a 
u 




1 1 "^00 


1 01 00 


1 9S00 


1 9^00 


1 0400 




1 0400 


1 0000 
xuuuu 








1 1 fiOO 


1 0/LOO 


1 4400 


1 9400 


1 1 fiOO 


1 0900 


SQOO 




2002cr 


-5 


9400 


10700 


9500 


12200 


11200 






9600 


9400 


2005am 


-4 




13200 


12000 


13400 


12900 


12100 




12100 


10600 


2001en 


-4 




14800 


13600 


13300 


13500 


12800 




11400 


10500 


2001br 


-2 




15700 


14500 


11900 


13700 


14100 




14200 


9100 


2002hw 


-1 




11900 


10700 


12900 


11900 


10200 




12200 




2005am 


+0 




13100 


11900 


11000 


13100 


10400 


11700 


12500 


7700 


2000dm 


+0 




11600 


10400 


9500 


11500 




10300 


10900 


6700 


2003du 


+ 1 




15000 


13800 


13500 


14900 


10800 


13700 


14200 




2001bf 


+4 




9800 


8600 


9700 


9600 


7400 


8500 


7600 


6800 


2000do 


+4 




11700 


10700 




11600 


7300 


9400 


10400 


3700 


2000dk 


+5 




10800 


9600 


9000 


10800 




9500 


11600 


2500 


2002ha 


+8 


10200 


11300 










10000 


11300 


3100 


2001bg 


+10 




12300 








10500 


11000 


12300 


3100 


2002ef 


+11 


9800 














11400 


2300 


2005am 


+14 




10500 










9200 


11400 


2900 


2004ab 


+18 




10300 








8200 


8900 


12200 


2600 


2002en 


+18 


12300 










7500 




13000 





Note. — 



These data are plotted in Figure [T9l 



-51 - 



Table 5. Estimated Line Strengths at 5,000K in the Region 0.7 — 2.4 /xm 



Ion 


Arest(air) 


No. lines 


log{ff/) 


Excitation 


Line Strength 


Line Strength 




fim 


in blend ^ 




Min. (eV) 


5,000K b 


10,000K b 




He i 


1.0830 


3 


-0.047 


19.82096 


15.1 


6.0 


He I 


2.0581 


1 


-0.424 


20.61718 


1.000 


1.000 


He I 


0.7065 


3 


-0.460 


20.96551 


0.410 


0.614 


ci 


1.0693 


7 


0.348 


7.48830 


1.000 


1.000 


Cl 


0.9093 


2 


0.142 


7.48830 


0.622 


0.622 


cl 


0.9406 


1 


0.225 


7.68529 


0.477 


0.599 


Cl 


0.9658 


1 


-0.264 


7.48830 


0.244 


0.244 


Cl 


1.4543 


1 


-0.110 


7.68529 


0.221 


0.277 


cl 


1.1754 


2 


0.661 


8.64774 


0.139 


0.535 


cl 


0.8335 


1 


-0.420 


7.68529 


0.108 


0.136 


cl 


1.1330 


1 


0.280 


8.53767 


0.075 


0.253 


cil 


0.6580 


2 


0.118 


14.44980 


88,393 


253 


ClI 


0.7235 


3 


0.330 


16.33422 


1,815 


46.3 


ClI 


1.8905 


2 


0.258 


19.49586 


1.000 


0.959 


ClI 


1.7846 


3 


0.550 


20.15184 


0.427 


0.878 


cll 


0.9903 


3 


1.010 


20.95206 


0.192 


1.000 


ClI 


0.9230 


3 


0.620 


20.84618 


0.100 


0.461 


ol 


0.7773 


3 


0.324 


9.14671 


1.000 


1.000 


ol 


0.8446 


3 


0.170 


9.52201 


0.294 


0.454 


ol 


0.9264 


9 


0.690 


10.74166 


0.057 


0.365 


ol 


1.1290 


9 


0.500 


10.98960 


0.021 


0.177 


ol 


0.7990 


6 


0.280 


10.98960 


0.013 


0.106 


ol 


1.3165 


3 


-0.095 


10.98960 


0.005 


0.045 


oil 


0.6695 


2 


-0.594 


23.44307 


2,155 


97.6 


oil 


2.1085 


2 


-1.690 


25.66290 


1.000 


0.595 


oil 


1.3811 


3 


-1.730 


25.66290 


0.912 


0.543 


oil 


1.4008 


2 


0.530 


29.62095 


0.017 


1.000 


oil 


0.7114 


3 


-0.170 


29.07061 


0.012 


0.378 


oil 


1.3029 


2 


0.460 


29.82223 


0.009 


0.674 


oil 


1.2500 


1 


0.210 


29.82223 


0.005 


0.379 


oil 


1.1667 


2 


0.870 


30.50576 


0.005 


0.784 


Mgl 


0.8807 


1 


-0.137 


4.34610 


1.000 


0.807 


Mgl 


1.1828 


1 


-0.290 


4.34610 


0.703 


0.568 


Mgl 


1.5033 


3 


0.340 


5.10817 


0.511 


1.000 


Mgl 


1.7109 


1 


0.140 


5.39409 


0.166 


0.453 


Mgl 


1.4878 


9 


0.660 


5.94632 


0.153 


0.790 


Mgl 


1.2083 


4 


0.450 


5.75363 


0.147 


0.609 


Mg ii 


0.9227 


2 


0.270 


8.65529 


1.000 


1.000 


Mg II 


1.0927 


3 


0.020 


8.86425 


0.346 


0.441 


Mg II 


0.7890 


3 


0.650 


10.00000 


0.106 


0.504 


Mg II 


0.8228 


2 


0.030 


10.00000 


0.025 


0.121 


Mg II 


1.0092 


3 


1.020 


11.63047 


0.006 


0.178 


Mg II 


0.9632 


3 


0.660 


11.56982 


0.003 


0.083 


Mg II 


2.1369 


1 


0.390 


11.50533 


0.002 


0.048 


Mg II 


1.6787 


3 


0.730 


12.08542 


0.001 


0.054 
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Table 5 — Continued 



Ion 


Arest(air) 


No. lines 


log(g/) 


Excitation 


Line Strength 


Line Strength 






in blend * 




Min. (eV) 


5,000K ^ 


10,000K ^ 




Sii 


1.2032 


1 


0.440 


4.95413 


1.000 


1.000 


Si I 


1.0827 


1 


0.220 


4.95413 


0.603 


0.603 


Si I 


1.0585 


1 


-0.020 


4.95413 


0.347 


0.347 


Si I 


1.5888 


1 


-0.030 


5.08269 


0.251 


0.292 


Si 11 


0.6355 


2 


0.297 


8.12157 


10,126 


49.6 


Si II 


1.6930 


2 


0.350 


12.14781 


1.000 


0.524 


Si II 


0.9413 


3 


0.980 


12.84018 


0.855 


1.000 


Si II 


0.7849 


3 


0.490 


12.52627 


0.573 


0.466 


Si II 


1.1737 


3 


0.620 


12.88100 


0.340 


0.416 


Si II 


1.3681 


2 


0.130 


12.88100 


0.110 


0.135 


Si II 


1.7183 


3 


0.890 


14.10541 


0.037 


0.187 


Si II 


1.1311 


3 


-0.460 


12.84018 


0.031 


0.036 


Si 111 


0.9324 


1 


-0.120 


20.55363 


7,223 


41.6 


Si III 


1.0526 


1 


-2.950 


20.55363 


10.7 


0.062 


Si III 


1.2523 


3 


0.500 


24.99684 


1.000 


1.000 


Si III 


0.7465 


6 


0.360 


24.99684 


0.724 


0.724 


Si III 


1.2601 


3 


0.160 


24.99662 


0.457 


0.457 


Si III 


0.9800 


1 


0.230 


25.33550 


0.245 


0.363 


Si III 


1.1341 


3 


0.430 


25.56430 


0.228 


0.441 


si 


0.9223 


3 


0.420 


6.52494 


1.000 


1.000 


si 


1.0457 


3 


0.260 


6.86061 


0.317 


0.469 


Si 


1.8940 


6 


0.450 


8.04620 


0.031 


0.183 


si 


2.2694 


3 


0.240 


7.87043 


0.029 


0.139 


si 


1.3809 


3 


0.110 


7.87043 


0.022 


0.103 


sll 


0.8315 


1 


-0.470 


14.06839 


1.000 


1.000 


SlI 


0.7967 


1 


-0.820 


14.00342 


0.519 


0.482 


SlI 


0.7590 


1 


-0.880 


14.23489 


0.264 


0.321 


SlI 


0.8855 


2 


-1.110 


14.16002 


0.185 


0.206 


SlI 


1.4501 


2 


-0.216 


16.59204 


0.005 


0.096 


SlI 


1.3529 


1 


-0.425 


16.53575 


0.004 


0.063 


Cal 


0.6149 


3 


0.100 


1.89906 


5.1 


2.0 


Cal 


0.7148 


1 


0.208 


2.70919 


1.000 


1.000 


Cal 


1.9753 


4 


-0.831 


1.89906 


0.599 


0.234 


Cal 


0.7326 


1 


0.073 


2.93271 


0.436 


0.565 


Cal 


1.9453 


1 


-1.094 


1.88593 


0.337 


0.130 


Call 


0.8538 


2 


-0.362 


1.70005 


13,944 


40.3 


Call 


0.8662 


1 


-0.623 


1.69252 


7,780 


22.3 


Call 


1.1839 


1 


0.300 


6.46831 


1.000 


0.731 


Call 


0.8921 


3 


0.729 


7.05003 


0.696 


1.000 


Call 


1.1950 


1 


0.000 


6.46831 


0.501 


0.367 


Call 


0.8249 


2 


0.621 


7.51535 


0.184 


0.454 


Call 


0.8202 


1 


0.315 


7.50564 


0.093 


0.227 


Call 


0.9931 


1 


0.072 


7.51535 


0.052 


0.128 


Mill 


1.2900 


1 


-1.059 


2.11436 


1.000 


0.541 


Mnl 


1.3310 


2 


-1.360 


2.14284 


0.468 


0.262 
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Table 5 — Continued 



Ion 


Ares* (air) 


No. lines 


log(3/) 


Excitation 


Line Strength 


Line Strength 




/im 


in blend ^ 




Min. (eV) 


5,000K b 


10,000K ^ 


Mnl 


1.3630 


3 


-1.516 


2.16386 


0.311 


0.178 


Mnl 


1.3859 


2 


-1.638 


2.17836 


0.227 


0.133 


Mnl 


1.5184 


2 


0.606 


4.88919 


0.074 


1.000 


Mnl 


0.8740 


3 


-0.055 


4.43521 


0.046 


0.370 


Mnl 


1.5263 


1 


0.379 


4.88919 


0.044 


0.593 


Mnl 


0.7309 


3 


-0.084 


4.43521 


0.043 


0.346 


Mnii 


0.7414 


3 


-2.202 


3.70608 


1.000 


1.000 


Mn II 


0.7342 


2 


-2.713 


3.70979 


0.306 


0.307 


Mn II 


0.9447 


1 


-2.389 


4.06545 


0.282 


0.428 


Mnll 


0.9388 


2 


-2.554 


4.07405 


0.189 


0.290 


Mnll 


0.8695 


1 


0.577 


9.24489 


0.002 


0.971 


Mnll 


0.8769 


1 


0.491 


9.24489 


0.001 


0.797 


Mn II 


0.8820 


1 


0.347 


9.24489 


0.001 


0.572 


Mn II 


0.9486 


2 


0.491 


9.46988 


0.001 


0.614 


Fc i 


0.8677 


2 


-1.212 


2.17609 


1.000 


0.417 


Fc I 


0.8824 


1 


-1.364 


2.19801 


0.670 


0.286 


Fc I 


1.1973 


1 


-1.476 


2.17609 


0.545 


0.227 


Fe I 


0.8388 


1 


-1.493 


2.17609 


0.524 


0.218 


Fel 


0.7511 


5 


0.107 


4.17798 


0.200 


0.851 


Fel 


0.8220 


2 


0.249 


4.32039 


0.199 


1.000 


Fel 


0.7196 


3 


-0.120 


4.10365 


0.141 


0.550 


Fc I 


0.7941 


2 


0.154 


4.38676 


0.137 


0.744 


Fe I 


0.7495 


2 


-0.102 


4.22065 


0.112 


0.500 


Fel 


0.7999 


1 


0.048 


4.37164 


0.111 


0.593 


Fell 


0.7712 


1 


-2.543 


3.90368 


1.000 


1.000 


Fe II 


0.7462 


1 


-2.734 


3.89187 


0.662 


0.653 


Fell 


0.9998 


1 


-1.826 


5.48450 


0.133 


0.832 


Fell 


1.0500 


2 


-1.997 


5.54914 


0.077 


0.521 


Fell 


1.0863 


1 


-2.121 


5.58957 


0.053 


0.374 


Fell 


1.1126 


1 


-2.236 


5.61560 


0.038 


0.278 


Fell 


0.7515 


1 


-2.362 


5.82361 


0.018 


0.163 


Fell 


0.9550 


1 


-2.027 


6.21915 


0.015 


0.223 


Felll 


0.7261 


1 


-2.910 


13.13443 


1.000 


0.3.54 


Felll 


0.7265 


1 


-1.708 


14.62496 


0.501 


1.000 


Felll 


0.9124 


1 


-2.577 


14.62496 


0.068 


0.135 


Felll 


1.6722 


1 


0.235 


22.87010 


0.000 


0.006 


FcIII 


0.7456 


1 


-1.905 


20.88295 


0.000 


0.000 


FcIII 


1.6672 


1 


0.057 


22.86637 


0.000 


0.004 


FcIII 


1.2039 


2 


0.663 


23.61137 


0.000 


0.007 


Felll 


1.2786 


1 


0.723 


23.67243 


0.000 


0.007 


Coi 


0.7085 


1 


-1.018 


1.88271 


1.000 


0.589 


Col 


0.7053 


1 


-1.440 


1.95586 


0.319 


0.205 


Col 


0.8098 


4 


0.290 


4.02115 


0.142 


1.000 


Col 


0.8027 


6 


0.116 


4.14628 


0.071 


0.579 


Col 


0.8374 


2 


-0.040 


4.07216 


0.059 


0.441 
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Table 5 — Continued 



Ion 


Arest(air) 


No. lines 


log{9/) 


Excitation 


Line Strength 


Line Strength 




/im 


in blend ^ 




Min. (eV) 


5,000K b 


10,000K ^ 




Co ii 


1.7772 


1 


-2.087 


5.04595 


1.000 


0.736 


Co II 


1.7462 


1 


-2.284 


5.12265 


0.532 


0.428 


Co II 


2.2205 


1 


-2.374 


5.04595 


0.516 


0.380 


Co II 


2.4596 


1 


-2.416 


5.12265 


0.392 


0.316 


Co II 


2.3613 


1 


-2.491 


5.17524 


0.292 


0.250 


Co II 


1.6064 


1 


-2.568 


5.12265 


0.276 


0.223 


Co II 


1.5759 


1 


-2.725 


5.04595 


0.230 


0.169 


Co II 


1.6361 


1 


-2.625 


5.17524 


0.215 


0.184 


Co II 


1.7239 


1 


-2.649 


5.17524 


0.203 


0.174 


Co II 


2.2497 


1 


-2.697 


5.20879 


0.168 


0.150 


Co II 


2.1347 


1 


-2.906 


5.04595 


0.152 


0.112 


Co II 


0.8342 


2 


0.753 


10.41714 


0.003 


1.000 


Co II 


0.8737 


1 


0.663 


10.41714 


0.002 


0.813 



^If more than one line is included in the blend, the indicated wavelength is a mean value, 
weighted by oscillator strength. Some lines have been ommitted from the table if their rest 
wavelengths are very close to that of a stronger line. 

The estimated line strength is computed using Nn f„ * 10"'^"/'=^, where /„ is the 
oscillator strength, g„ is the statistical weight, Xn is the excitation potential, and T is the 
temperature. The values are normalized by dividing by the value for one of the strongest 
lines for each ion. This is usually the strongest line at that temperature, but in some cases 
(for example Ca II) the strongest line is so much stronger that the normalization is nor 
relevant. In that case we chose the second or third strongest line, (see i|3.1ll 
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Table 6. Estimated Line Strengths at 10,000K in the Region 0.7 - 2.4 /xm 



Ion 


Arest(air) 


No. lines 


log{ff/) 


Excitation 


Line Strength 


Line Strength 




fim 


in blend ^ 




Min. (eV) 


5,000K b 


10,000K b 




He i 


1.0830 


3 


-0.047 


19.82096 


15.1 


6.0 


He I 


2.0581 


1 


-0.424 


20.61718 


1.000 


1.000 


He I 


0.7065 


3 


-0.460 


20.96551 


0.410 


0.614 


ci 


1.0693 


7 


0.348 


7.48830 


1.000 


1.000 


Cl 


0.9093 


2 


0.142 


7.48830 


0.622 


0.622 


cl 


0.9406 


1 


0.225 


7.68529 


0.477 


0.599 


Cl 


1.1754 


2 


0.661 


8.64774 


0.139 


0.535 


Cl 


1.7325 


8 


1.030 


9.70241 


0.028 


0.368 


cl 


1.6890 


1 


0.568 


9.00319 


0.049 


0.286 


cl 


1.4543 


1 


-0.110 


7.68529 


0.221 


0.277 


cl 


1.1330 


1 


0.280 


8.53767 


0.075 


0.253 


cil 


0.6580 


2 


0.118 


14.44980 


88,393 


253 


cll 


0.7235 


3 


0.330 


16.33422 


1,815 


46.3 


ClI 


0.9903 


3 


1.010 


20.95206 


0.192 


1.000 


ClI 


1.8905 


2 


0.258 


19.49586 


1.000 


0.959 


cll 


1.7846 


3 


0.550 


20.15184 


0.427 


0.878 


ClI 


0.9230 


3 


0.620 


20.84618 


0.100 


0.461 


ol 


0.7773 


3 


0.324 


9.14671 


1.000 


1.000 


ol 


0.8446 


3 


0.170 


9.52201 


0.294 


0.454 


ol 


0.9264 


9 


0.690 


10.74166 


0.057 


0.365 


ol 


1.1290 


9 


0.500 


10.98960 


0.021 


0.177 


ol 


1.8021 


18 


0.880 


12.07943 


0.004 


0.120 


ol 


0.7990 


6 


0.280 


10.98960 


0.013 


0.106 


Mg 11 


0.9227 


2 


0.270 


8.65529 


1.000 


1.000 


Mg II 


0.7890 


3 


0.650 


10.00000 


0.106 


0.504 


Mg II 


1.0927 


3 


0.020 


8.86425 


0.346 


0.441 


Mg II 


1.0092 


3 


1.020 


1 1 63047 


0.006 


0.178 


Mg II 


0.8228 


2 


0.030 


10.00000 


0.025 


0.121 


Mg II 


0.9632 


3 


0.660 


11.56982 


0.003 


0.083 


Mg II 


1.8606 


9 


1.220 


12.85873 


0.001 


0.068 


Mg II 


1.6787 


3 


0.730 


12.08542 


0.001 


0.054 


Mg II 


2.1369 


1 


0.390 


11.50533 


0.002 


0.048 


Si 11 


0.6355 


2 


0.297 


8.12157 


10,126 


49.6 


Si II 


0.9413 


3 


0.980 


12.84018 


1.000 


1.000 


Si II 


1.6930 


2 


0.350 


12.14781 


1.169 


0.524 


Si II 


0.7849 


3 


0.490 


12.52627 


0.671 


0.466 


Si II 


1.1737 


3 


0.620 


12.88100 


0.397 


0.416 


Si II 


1.7183 


3 


0.890 


14.10541 


0.043 


0.187 


Si II 


1.3681 


2 


0.130 


12.88100 


0.128 


0.135 


Si II 


2.1967 


3 


0.680 


14.13227 


0.025 


0.112 


Si 111 


0.9324 


1 


-0.120 


20.55363 


7,223 


41.6 


Si III 


1.2523 


3 


0.500 


24.99684 


0.094 


1.000 


Si III 


0.7465 


6 


0.360 


24.99684 


0.068 


0.724 


Si III 


1.2601 


3 


0.160 


24.99662 


0.043 


0.457 


Si III 


1.1341 


3 


0.430 


25.56430 


0.021 


0.441 



-56- 



Table 6 — Continued 



Ion 


Ares* (air) 


No. lines 


log(ff/) 


Excitation 


Line Strength 


Line Strength 




fim 


in blend ^ 




Min. (eV) 


5,000K b 


lO.OOOK ^ 


Si III 


0.9800 


1 


0.230 


25.33550 


0.023 


0.363 


Si III 


0.8266 


6 


0.830 


26.65707 


0.004 


0.311 


Si III 


0.7612 


1 


0.660 


26.60039 


0.003 


0.225 


Si III 


1.5015 


1 


0.030 


25.77484 


0.005 


0.137 


Si III 


0.9375 


1 


-0.490 


25.33550 


0.004 


0.069 


Si III 


0.7729 


2 


-0.680 


24.99669 


0.006 


0.066 


Si III 


1.0526 


1 


-2.950 


20.55363 


1.000 


0.062 


Sii 


0.8315 


1 


-0.470 


14.06839 


1.000 


1.000 


SlI 


0.7967 


1 


-0.820 


14.00342 


0.519 


0.482 


SlI 


0.7590 


1 


-0.880 


14.23489 


0.264 


0.321 


sll 


0.7721 


1 


-0.990 


14.29456 


0.179 


0.232 


sll 


0.8855 


2 


-1.110 


14.16002 


0.185 


0.206 


sll 


1.4501 


2 


-0.216 


16.59204 


0.005 


0.096 


sll 


1.3529 


1 


-0.425 


16.53575 


0.004 


0.063 


Call 


0.8538 


2 


-0.362 


1.70005 


13,944 


40.3 


Call 


0.8662 


1 


-0.623 


1.69252 


7,780 


22.3 


Call 


0.8921 


3 


0.729 


7.05003 


0.696 


1.000 


Call 


1.1839 


1 


0.300 


6.46831 


1.000 


0.731 


Call 


0.8249 


2 


0.621 


7.51535 


0.184 


0.454 


Call 


1.1950 


1 


0.000 


6.46831 


0.501 


0.367 


Call 


0.9891 


3 


0.943 


8.43855 


0.045 


0.327 


Call 


0.8202 


1 


0.315 


7.50564 


0.093 


0.227 


Call 


1.8850 


2 


0.878 


9.01809 


0.010 


0.144 


Call 


0.9931 


1 


0.072 


7.51535 


0.052 


0.128 


Mn 11 


0.7414 


3 


-2.202 


3.70608 


1.000 


1.000 


Mn II 


0.8695 


1 


0.577 


9.24489 


0.002 


0.971 


Mnll 


0.8769 


1 


0.491 


9.24489 


0.001 


0.797 


Mnll 


0.9486 


2 


0.491 


9.46988 


0.001 


0.614 


Mnll 


0.8820 


1 


0.347 


9.24489 


0.001 


0.572 


Mn II 


0.9447 


1 


-2.389 


4.06545 


0.282 


0.428 


Mn II 


0.7342 


2 


-2.713 


3.70979 


0.306 


0.307 


Mnll 


0.9388 


2 


-2.554 


4.07405 


0.189 


0.290 


Mnll 


0.7221 


3 


0.768 


10.67048 


0.000 


0.288 


Mnll 


0.9444 


1 


0.129 


9.46988 


0.000 


0.267 


Mnll 


0.8110 


1 


-2.002 


■T.. 37762 


0.033 


0.228 


Mn II 


1..5101 




o.2:-!7 


<).8()()00 


0.000 


0.216 


Fell 


0.8288 


1 


0.693 


9.65426 


0.003 


1.000 


Fell 


0.9297 


1 


0.412 


9.65426 


0.001 


0.524 


Fell 


0.7712 


1 


-2.543 


3.90368 


1.000 


0.460 


Fell 


0.9998 


1 


-1.826 


5.48450 


0.133 


0.382 


Fell 


0.9095 


1 


0.269 


9.65426 


0.001 


0.377 


Fell 


0.9183 


2 


0.245 


9.70089 


0.001 


0.338 


Fell 


0.7462 


1 


-2.734 


3.89187 


0.662 


0.300 


Fell 


1.0500 


2 


-1.997 


5.54914 


0.077 


0.239 


Fell 


1.0863 


1 


-2.121 


5.58957 


0.053 


0.172 
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Table 6 — Continued 



Ion 


Arest(air) 


No. lines 


log{3/) 


Excitation 


Line Strength 


Line Strength 




/im 


in blend ^ 




Min. (eV) 


5,000K ^ 


10,000K ^ 




Fein 


0.7265 


1 


-1.708 


14.62496 


0.501 


1.000 


Felll 


0.7261 


1 


-2.910 


13.13443 


1.000 


0.354 


Felll 


0.9124 


1 


-2.577 


14.62496 


0.068 


0.135 


Felll 


1.2786 


1 


0.723 


23.67243 


0.000 


0.007 


Felll 


1.2039 


2 


0.663 


23.61137 


0.000 


0.007 


Felll 


1.6722 


1 


0.235 


22.87010 


0.000 


0.006 


Felll 


1.2955 


2 


0.602 


23.67243 


0.000 


0.006 


Co 11 


0.8342 


2 


0.753 


10.41714 


0.003 


1.000 


Co 11 


0.8737 


1 


0.663 


10.41714 


0.002 


0.813 


Co 11 


1.7772 


1 


-2.087 


5.04595 


1.000 


0.736 


Co 11 


0.8582 


1 


0.598 


10.48815 


0.002 


0.644 


Co 11 


0.9070 


1 


0.433 


10.41714 


0.001 


0.479 


Co 11 


1.7462 


1 


-2.284 


5.12265 


0.532 


0.428 


Co 11 


0.8803 


2 


0.521 


10.70637 


0.001 


0.419 


Co 11 


2.2205 


1 


-2.374 


5.04595 


0.516 


0.380 


Co 11 


0.8512 


2 


0.426 


10.61330 


0.001 


0.375 


Co 11 


2.4596 


1 


-2.416 


5.12265 


0.392 


0.316 


Co 11 


2.3613 


1 


-2.491 


5.17524 


0.292 


0.250 



^If more than one line is included in the blend, the indicated wavelength is a mean value, 
weighted by oscillator strength. Some lines have been ommitted from the table if their rest 
wavelengths are very close to that of a stronger line. 

°The estimated line strength is computed using Nn U * lO-X"/''^, where /„ is the 
oscillator strength, g-n is the statistical weight, Xn is the excitation potential, and T is the 
temperature. The values are normalized by dividing by the value for one of the strongest lines 
for each ion. This is usually the strongest line at that temperature, but in some cases (for 
example Ca 11) the strongest line is so much stronger that the normalization is nor relevant. 
In that case we chose the second or third strongest line, (see i|3.1l l 
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Fig. 1. — The spectrum from SN 2000dm at +0 day provides an example of our spike 
elimination routine (see text The original data is plotted in black and the cleaned 

spectrum is superimposed in blue (see discussion 
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5 10 15 20 25 30 35 40 45 50 10 20 30 40 50 

Frequency Bin Frequency Bin 

Fig. 2. — The power spectrum and filter used in our reduction package (see text §2.11) . The 
left panel is a log plot of the first fifty points (out of 1901) in one half of the power spectrum 
for the wavelength region 0.80-1.35 fim from the SXD data from SN 2005am at -4 days. The 
noise level (N) is taken to be a constant equal to the mean of the next 100 points (50:150). 
The location of N is indicated in the figure by a horizontal dashed line. The signal (S) 
descends from near -4 on the left edge and intersects the noise level at point 27. The right 
panel shows the filter produced by the formula: S/(S+N) (see discussion §2.ip . 
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Frequency Bin 



1500 



2000 



Fig. 3. — Tlie Fourier transform for the wavelength region 0.80-1.35 fim from SN 2005am 
at -4 days (see text §2.ip . The initial frequency spectrum is plotted in black. The same 
spectrum after processing with the filter as described in the text is superimposed in red (see 
discussion §2.11) . 
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0.4 r 




— 0.8 H , , ^ \ ^ ^ ^ \ ^ ^ , I , , ^ \ ^ ^ ^ \ ^ , , 

0.84 0.86 0.88 0.90 0.92 0.94 0.96 

Rest Wavelength of Host Galaxy (/xm) 




— 0.3 E , , , \ i i i \ i i , \ , , , I i i i \ i i i I , , , z 

0.98 1.00 1.02 1.04 1.06 1.08 1.10 1.12 
Rest Wavelength of Host Galaxy (/xm) 

Fig. 4. — Two absorption features from tfie liigli S/N SXD spectrum of SN 2005am at -4d 
demonstrate that tlie FT smoothing technique is robust and the measured velocities do not 
change significantly due to fine interpretations of the smoothing parameter by the user (see 
text §2.11 and Table [1]). The raw data are black and the smoothed spectrum in black is the 
Fourier transform spectrum as with a smoothing parameter of 27 as shown in Figure [2^. 
The spectra produced by other smoothing parameters are plotted in different colors (see 
discussion §2.ip . On this scale, the distance between ticks on the wavelength axis (0.005 
yum) is about 1,500 km s~^. 
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0.8 0.9 1.0 1.1 1.2 1.3 

Rest Wavelength of Host Galaxy (fim) 



Fig. 5. — Two spectra with different S/N show how use of the FT smoothed spectra facihtate 
understanding of the overall spectral shape and identification of the precise locations of 
maxima and minima (see text §2.ip . The raw data are plotted in green and the FT smoothed 
spectra are superimposed in black (see discussion §2.ip . 
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0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.40.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Rest Wavelength of Host Galaxy (/xm) 



Fig. 6. — Raw data from four spectra showing possible fits to the continua (see text §2.2p . 
The top left spectrum is LRS, the others are SXD. The red line is a -3.0 power law and the 
green line is a -4.5 power law. It is clear that a single power law does not fit any of the 
spectra (see discussion §2.2p . 
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I \ \ \ \ \ \ \ I \ \ \ \ I \ \ \ \ I \ \ \ I \ \ \ \ I \ \ \ \ I \ \ \ \ \ \ \ \ I \ \ \ \ I \ \ \ \ \ \ \ \ I \ \ 

0.8 0.9 1.0 1.1 1.2 1.3 

Rest Wavelength of Host Galaxy (/xm) 

Fig. 7. — Estimated fits for local continua in the vicinity of absorption features are plotted 
on the figure in red, green, and blue (see text §2.2p . The top spectrum is a raw data SXD 
spectrum from SN 2001en at -4 days. On the bottom is a raw LRS spectrum from SN 2002cr 
at -5 days. These are the same spectra displayed in the top two panels of Figure M and the 
fitted regions are expanded in Figure M (see discussion §2.21) . 
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1.00 1.02 1.04 1.06 1.06 1.10 1.00 1.02 1.0-1 1.06 l.OB 1.10 




1.16 1.18 1.20 1.22 1.24 1.26 1.16 I.IB 1.20 1.22 1.24 1.26 

Rest Wavelength of Host Galaxy {fim) 

Fig. 8. — Expanded spectral regions from Figure [7] with the continuum normahzed to zero 
(see text §2.2p . The left column displays SXD data from SN 2001en (-4d) and the right 
column is LRS data from SN 2002cr (-5d). The top row corresponds to the regions in Figure 
[7] marked with a green continuum, the middle row is from the red regions and the bottom 
row is from the blue regions. The FT smoothed spectrum is superimposed on the raw data 
in red (see discussion §2.2p . 
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Fig. 9. — Forty-one spectra from normal SNe la obtained between -14d and +75d with 
respect to V^ax (see text The abscissa is wavelength in microns, shifted to the rest 
frame of the host galaxy and the ordinate is log flux. All spectra have been normalized to 
1.0 at 1.0 fim and shifted by a constant for clarity. Spectra plotted in yellow have low S/N 
and will be omitted for the analysis of spectral features. Observational and spectral details 
are provided in Tables [2] and [3l 
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1.0 1-2 1.4 1.6 l.B 2.0 2.2 2.4 
Rest Wavelength of Host Galaxy (ptrn) 



Fig. 10. — Thirty-four spectra selected from our sample for S/N ratio greater than 10 in the 
region 1.10 — 1.30 /im (see Tables [2] and [3]) . These spectra cover the entire temporal range 
of our sample, from -14 to -|-75 days with respect to V^ax (see text The abscissa is 
wavelength, in microns, shifted to the rest frame of the host galaxy and the ordinate is log 
flux. All spectra have been normalized to 1.0 at 1.0 /im and shifted by a constant for clarity. 
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^ I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Rest Wavelength of Host Galaxy (/xrn) 

Fig. 11. — A very early spectrum of SN la 2002fk obtained fourteen days before Vmax', only 
5 or 6 days after the explosion. The data are displayed in both linear and log space (plus 
a constant) to demonstrate how the sizes of features at longer wavelengths are exaggerated 
by the log scale. Features are labeled for which we have reasonably confident identifications. 
A general discussion of this spectrum can be found in §3.1.11 and a more detailed analysis is 
found in lA.li 
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Fig. 12. — Eighteen spectra from the Photospheric Phase of SNe la obtained between -14d 
and +5d relative to Vmax- During this epoch the photosphere is expected to be in the outer 
layers of the SN and the line-forming regions will be close to the photosphere. The figure 
displays only the wavelength region 0.65 - 1.35 fim for detail. The complete spectra are 
found in Figure [TOl A general discussion of features detected during the Photospheric Phase 
can be found in §3.1.21 and a more detailed analysis in found in IA.2I 
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Fig. 13. — Four representative spectra from the Photospheric Phase of spectral development 
with features labeled. The spectra are from SNe 2004bw (LRS obtained at -9d), 2005am 
(SXD obtained at -4d), 2000dm (SXD obtained at +0d), and 2001bf (SXD obtained at +4d). 
The top panel displays only the wavelength region 0.65 - 1.35 /im for detail while the bottom 
panel contains the complete spectra. A general discussion of these features can be found in 
§3.1.21 and a more detailed analysis in found in IA.2[ 
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0.6 0,8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 
Rest Wavelength of Host Galaxy (/^m) 



Fig. 14. — Ten spectra from the Extended Photospheric Phase of SNe la obtained between 
+4d and +24d with respect to to Vmax- During the phase represented by this group, the shape 
of the continuum becomes influenced by hne-blanketing from Fe-group hues which extends 
the effective photospheric radius (see §5.3p . A general discussion of features detected during 
the Extended Photospheric Phase can be found in §3.1.31 and a more detailed analysis is 
found in §A.3[ 
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0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 

Rest Wavelength of Host Galaxy (/^m) 

Fig. 15. — Four representative spectra from the the Extended Photospheric Phase of spectral 
development with features labeled. The spectra are from SNe 2000dk (+5d), 2001bg (+10d), 
2004ab (+18d), and 2004da (+24d). The top panel displays only the wavelength region 0.65 
- 1.35 fim for detail while the bottom panel contains the complete spectra. A general 
discussion of these features can be found in §3.1.31 and a more detailed analysis in found in 
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Rest Wavelength of Host Galaxy (jOtm) 

Fig. 16. — Seven spectra from the Transitional Phase SNe la obtained between +37d and 
+75d with respect to to Vmax- This epoch is a transitional phase between the photospheric 
era and a true nebular phase. This group displays a lot of new structure relative to earlier 
NIR spectra and the features are remarkably similar from one spectrum to the next. A 
general discussion of features detected during the Co II Phase can be found in §3.1.41 and a 
more detailed analysis is found in §A.4[ 
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Fig. 17. — Fifteen spectra between +14d and +75d are presented with the positions of strong 
Co II and Fe II hnes drawn on the figure so that the indicated velocities are 11,000 km s~^ 
for the top spectrum and 6.000 km s~^ for the bottom spectrum. This corresponds to the 
measured velocities for features from spectr a during this epoc h. The data plotted in yellow 
(SNe 1998aq and 1998bu) were obtained by iRudy fc Puetterl and have not been previously 
published. Co II lines are easily detected in the all spectra in this group and and Fe II lines 
begin to appear about one month post-maximum (see §3.1.31 and 13.1.41) . 
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Fig. 18. — The uniformity of feature evolution in NIR spectra of SNe la within a few weeks 
of maximum hght is demonstrated by interleaving "snapshot" data from 19 SNe la in our 
sample with two time sequences of spectra obtained from individual SNe la (see text §4.ip . 
Spectral shape and feature development with time are found to be consistent in individual 
events and also between diver se SNe la. Five spec tra from SN 1999ee obtained between -lid 
and +13d are shown in red (IHamuy. et al.ll2002l ). The four spectra in blue are from our 
sample and show SN 2005am between -4d and +14d. 
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Fig. 19. — Plot of measured Doppler velocities for ions with confident identifications in 
spectra obtained before +20d. In all panels, the ordinate is the time of observation in days 
relative to Vmax, and the abscissa is the Doppler velocity in 10^ km s^^. The O I data 
include measurements from two lines (0.7774 and 0.9264 /im), the Mg II data include three 
lines (0.9227, 1.0092, and 1.0927 /im), and there is one line each from Si II (0.9413 /im), 
S I (0.9223 /im), the Ca II IR-triplet (0.8538 /zm blended line), and Fe II (0.9998 /im). The 
measurements are compiled in Table H] and discussed in §4.21 
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Fig. 20. — Eighteen NIR spectra from SNe la as in Figure [I2] (-14d to +5d). The locations of 
strong C I hnes marked on the figure at Doppler velocities of 11,000 km s~^ which is typical 
for O I during this epoch. The strongest lines, as indicated in Table O are marked in red and 
the weaker lines appear in black. Small features appear in some spectra near the expected 
location for lines from 0.9406, 1.1330, and 1.1754 /im, but if these lines produce absorption 
features then the stronger lines at 1.0693 and 0.9093 fim should be detected and they are 
not. The possible feature for the 0.9406 fim line is more likely attributed to Si II at 0.9413 
/im (see §4.3p . 
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Fig. 21. — Model data for radius at r = 1.0 as a function of wavelength compared to spectra 
obtained from SN 2005am at epochs: -4d, +0d, +7d, and +14d. The models clearly predict 
the observed behavior of a pseudo photosphere created at extended radii that increases 
the continuum flux in wavelength regions: 0.9-1.2, 1.5-1.85, and 2.05-2.45 /im. Increased 
opacity due to line-blanketing from thousands of Fe-group lines increases the effective radius 
of the photosphere at certain wavelengths. The model resolution creates individual peaks 
but groups of lines will be blended by Doppler shifting to smear out the sharp features, (see 
text CT) . 
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